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ABSTRACT 


A  series  of  computer  programs  is  under  development 
for  use  in  the  design  and  evaluation  of  roll  stabiliza¬ 
tion  devices.  This  report  is  the  user  manual  for  the 
first  program  that  has  been  completed.  Identified  by 
the  acronym  FINCON,  the  program  is  based  on  the  work  of 
Cox  and  Lloyd  published  in  Volume  85  of  the  Transactions 
of  the  Society  of  Naval  Architects  and  Marine  Engineers. 
FINCON  predicts  stabilized  and  unstabilized  ship  roll 
motion,  bilge  keel  and  antiroll  fin  sizing  effects,  and 
the  influence  of  fin  controller  characteristics  by  use 
of  a  one  degree-of-freedom  roll-motion  equation.  Non¬ 
linear  roll  damping  characteristics,  derived  from  model 
experiments  or  by  other  means,  are  incorporated  by  a  com¬ 
bination  of  equivalent  linearization  and  an  iteration 
procedure.  Results  are  predicted  for  short-crested  seas 
(for  stabilizer  design  purposes)  and  long-crested  seas, 
which  are  described  by  two-parameter  Bretschnieder  wave 
spectra. 


ADMINISTRATIVE  INFORMATION 

The  development  and  documentation  of  the  computer  program  reported 
herein  is  a  part  of  the  Conventional  Ship  Seakeeping  Research  and  Develop¬ 
ment  Program  (Block  SF  43  421  202)  and  the  Ship  Performance  and  Hydrodynam¬ 
ics  Program  (Block  SF  43  421  001)  both  under  Program  Element  62543.  At 
David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  it  is 
identified  by  Work  Units  1504-100,  1507-200,  and  1500-104.  Authors  Susan 
L.  Bales  and  Geoffrey  G.  Cox  are  DTNSRDC  personnel.  Author  John  R.  Tucker 
is  on  the  staff  of  Chi  Associates,  Inc. 

INTRODUCTION 

Devices  such  as  bilge  keels,  anti-roll  fins,  and  anti-roll  tanks  have 
been  used  over  the  years  to  reduce  the  roll  motion  of  naval  and  commercial 
vessels.  In  recent  years,  the  U.S.  Navy  has  become  increasingly  involved 
in  the  design  and  development  of  suitable  roll  stabilization  devices  for 
Navy  ships.  This  report  provides  a  user's  manual  for  a  computer  program 
which  permits  prediction  of  unstabilized  and  bilge  keel/anti-roll,  fin- 
stabilized,  ship  roll  motions;  bilge-keel  and  fin-sizing  requirements;  and 
the  influence  of  fin  controller  characteristics.  The  program  is  known  by 
the  acronym  FINCON. 
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The  need  for  improved  design  and  performance  evaluation  tools  for 
bilge  keels,  antiroll  fins,  and  their  controllers  has  been  recognized  by 
Cox  and  Lloyd,  who  provide  the  hydrodynamic  basis  for  such  Investigations. 
A  comprehensive  work.  Reference  1  covers  such  topics  as  the  current  state- 
of-the-art  for  bilge  keels,  fins  and  tanks,  the  status  of  lateral  motion 
predictions,  measures  of  effectiveness,  design  practices,  and  sea  state 
specifications  for  design.  The  computer  program  described  in  this  report 
is  based  on  the  procedures  detailed  in  Reference  1. 

Reference  1  recognizes  the  potential  need  for  roll  and  roll  stabiliza¬ 
tion  prediction  tools  throughout  the  so-called  design  spiral  of  U.S.  Navy 
ships.  Hence,  the  procedures  outlined  there  and  used  here  can  be  employed 
with  a  very  simple  descriptor  of  ship  particulars  and  geometry.  Specifi¬ 
cally,  modifications  to  the  one-degree-of-f reedom  roll  motion  equation  of 
2 

Conolly  are  used  to  predict  ship  roll  motion,  and  the  required  program 
input  is  rather  easy  to  obtain  from  the  data  available  during  early  stages 
of  ship  design.  For  instance,  in  addition  to  specification  of  sea  condi¬ 
tion  and  ship  speed,  particulars  such  as  ship  length,  beam  displacement, 
transverse  metacentric  height  and  radius,  natural  roll  period,  and  roll 
decay  coefficients  are  required.  The  estimates  of  these  required  input 
variables  can  be  refined  as  the  design  process  continues  and  more  accurate 
data  becomes  available  through  model,  experiments,  etc.  The  manpower  and 
computer  time  costs  involved  in  such  predictions  are  relatively  low,  and 
these  predictions  can  be  completed  quickly  in  comparison  to  other  sea¬ 
keeping  design  evaluation  procedures. 

An  added  feature  of  the  FINCON  program  is  the  capability  to  recognize 
the  effect  of  fin  saturation,  which  occurs  in  heavy  seas,  on  RMS  (root  mean 
square)  roll  angle.  The  approach  is  based  on  a  refinement  of  the  method 
given  in  Appendix  3  of  Reference  1,  and  specific  details  of  this  improve¬ 
ment  will  be  published  in  a  future  report  currently  under  preparation  by 
Cox. 


*A  complete  listing  of  references  is  given  on  page  59. 
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PROGRAM  ORGANIZATION 

FINCON,  written  in  extended  FORTRAN,  is  operable  on  DTNSROC's  CDC 
6000  computers. 

Figure  1  is  an  overview  of  the  organization  of  FINCON.  The  overall 
program  consists  of  two  major  parts,  FINCON  and  FINSTAB  with  FINCON  acting 
as  the  driver  or  so-called  main  program. 


IESSJ  BRWSSP  ITREQ  I  ALGRNG  I  I  ERRF 


Figure  1  -  Program  Organization 

As  the  driver  of  the  rest  of  the  program,  FINCON  calls  all  other 
routines  for  execution.  In  addition,  FINCON  provides  the  unstabilized  roll 
calculations  (including  the  effect  of  bilge  keels)  and  is  responsible  for 
all  program  input  and  output. 

FINSTAB  is  the  second  major  part  of  the  program.  It  predicts  fin- 
stabilized  roll  angles,  fin  angles,  fin  angular  velocities,  and  fin  angular 
accelerations  for  the  input  fin  and  fin  controller  characteristics. 

FINSTAB  also  includes  the  fin  saturation  calculations.  Although  it  is 
celled  into  execution  by  FINCON,  given  the  proper  input  and  output  speci¬ 
fications,  FINSTAB  in  reality  can  be  executed  as  an  independent  program. 


BESSJ  Is  a  COC  system  subroutine  which  is  called  up  at  the  time  of 
execution  for  calculation  of  Bessel  functions  of  the  first  kind.  BESSJ  is 
only  used  when  the  ship's  waterplane  is  declared  elliptical. 

BRWSSP  is  a  subroutine  to  predict  the  Bretschnieder  wave-slope  spectrum 
for  a  specified  significant  wave  height  and  modal  wave  period.  Initially, 
the  wave-height  spectrum  is  computed;  and  then  it  is  converted  to  a  wave- 
slope  spectrum  by  multiplying  by  the  product  of  a  constant  and  the  square 
of  the  wave  number  k, 

k  =  w2/g  (1) 

where  is  the  wave  frequency  in  radians  per  second  and  g  is  the  accelera- 

2 

tion  due  to  gravity.  The  constant  (180/ir)  enters  in  to  permit  conversion 

into  degrees  to  yield  values  of  the  roll-response  amplitude  operator,  which 

2 

is  computed  by  FINCON  in  units  of  (degree/degree)  . 

ITREQ  is  a  subroutine  which  iterates  between  an  equivalently  linear¬ 
ized  roll  damping  curve,  which  is  a  function  of  RMS  roll  rate,  and  the 
predicted  RMS  roll  rate.  In  brief,  the  iteration  continues  until  the 
computed  short-  or  long-crested  RMS  roll  rate,  either  unstabilized  or 
stabilized,  is  within  a  small  value,  epsilon,  of  the  previously  computed 
value.  The  appropriate  RMS  roll  angle  can  then  be  found.  Additional 
details  of  the  exact  procedure  are  provided  in  the  Appendixes  A,  B,  and  C. 
ITREQ  is  called  by  both  FINCON  and  FINSTAB. 

ERRF  is  a  function  subprogram  which  gives  a  rational  approximation  to 
the  error  integral 


and  is  used  in  the  calculation  of  fin  saturation  effects. 

ALGRNG  is  an  integration  subroutine  which  performs  a  so-called 
Lagrangian  or  quadratic  integration  over  three  points  at  a  time.  The 
subroutine  is  called  by  both  FINCON  and  FINSTAB,  and  spectral  closure  is 
ensured  by  the  techniques  outlined  in  Figure  2,  adopted  from  Reference  3. 

•  4 


□  A, 


AREA  OF  COMPUTED  RESPONSE  SPECTRUM 
1/2S1  WHERE  wQ  -  w,  -  0.03 

AREA  OF  RIGHT  TRIANGLE  FORMED  BY 
DRAWING  STRAIGHT  LINE  THROUGH  $. 
AND  S2  TO  THE  ABSCISSA.  TRIANGLE 

SAME  AS  A3  BUT  FOR  S„.1  AND  S„ 


0.05 

2.00 


«  (OR  Wg) 


IF  S1  >0.1  S^,  THEN  A,  -A,  +  MINIMUM  OF  A2  AND  A3 
IF  Sn<Sn_1AND  S„>0.1  Smx,  THEN  A,<-A1+A4 


Figure  2  -  Spectral  Closure  Procedure 
(From  Reference  3) 
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The  procedure  for  predicting  the  unstabilized  and  stabilized  roll 
angles,  as  well  as  fin  angles,  and  velocities  is  outlined  in  Appendix  A. 
Because  the  procedure  is  developed  in  detail  in  Reference  1,  only  a  listing 
of  the  equations  used  is  given  in  Appendix  A.  In  general,  the  FORTRAN 
variables  have  been  named  as  closely  as  possible  to  the  variables  of  the 
equations  of  Appendix  A,  so  the  user  should  have  relatively  little  diffi¬ 
culty  in  "reading"  the  program.  A  listing  of  the  program  is  given  in 
Appendix  B.  Appendix  C  describes  two  special  algorithms  used,  namely  the 
iteration  algorithm  and  the  cosine  squared  short-crested  sea  algorithm  for 
arbitrary  spreading  angles. 


PROGRAM  INPUT 

As  described  in  the  Introduction,  the  input  requirements  of  F1NC0N 
are  relatively  simple.  The  input  variables  consist  of  sea  conditions,  ship 
speed,  and  very  simple  descriptions  of  the  ship  particulars  and  geometry. 
Unlike  many  multi-degree-of-freedom  ship-motion  prediction  programs  in  use 
today,  neither  offsets  of  the  ship  sections  nor  Lewis  forms  are  required. 
This  makes  implementation  of  FINCON  possible  at  an  early  design  stage 
(e.g.,  before  the  ship  lines  are  "firmed  up"). 

Table  1  describes  the  required  input  to  the  program  and  Table  2 
explains  the  notation  and  variable  names.  Up  to  eight  sea  conditions  and 
five  ship  speeds  may  be  executed  within  a  single  run,  using  either  metric 
or  English  units  for  input/output.  If  ITERATE  on  card  10  is  greater  than 
zero,  then  roll  damping  is  treated  as  nonlinear,  and  coefficients,  DUC, 
should  be  input.  For  example,  letting  d  *  DUC  for  a  simpler  notation,  then 

n  -  dQ  +  1.61dpyP  +  1.88^  +  4d2y2  +  9.4d3y3  +  24d4y4  (3) 

where  y  «*  O^/u)^  is  the  RMS  roll  rate  divided  by  ship  natural  frequency, 
p  is  0.772,  and  n  is  the  roll  damping  coefficient  as  a  function  of  y. 
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TABLE  1  -  INPUT  TO  FINCON  PROGRAM 


Contents 


NAME1,  NAME2,  NAME3 
TITLE2 

NWH,  NV,  I UNITS 

(SWH(IWH),  IWH-1,  NWH)  where 
1  £  NWH  £  10 

(TO(IWH),  IWH-1,  NWH)  where 
1  £  NWH  £  10 

(VK(IV),  IV-1,  NV)  where 
1  £  NV  £  5 

DISPTON,  L,  T,  GM,  BM,  TPHI,  Q 
SHAPE 

ISC,  ANGLE 
ITERATE 

IF  ITERATE  =  0 

((DUC(IV,1) ,  IV-1,  NV)  where 

1  <  NV  <  5 


IF  ITERATE  #  0 

((DUG (IV, 1) ,  1-1,  6),  IV-1,  NV) 
where  1  £  NV  £  5 

IPRINT(l) ,  IPRINT(2) 

NSTAB 


Format 


(315) 

(8F10.5) 

(8F10.5) 

(5F10.5) 

(7F10.5) 

(A10) 

(I5.F10.5) 

(15) 

(5F10.5) 


(6F10.5) 

(2A10) 

(15) 


Repeat  cards  14  through  17  NSTAB  times. 


M,  AREA,  R 

(I5.2F10.5) 

DCLDBFS ,  HO,  HI,  H2,  H3,  H4 

(6F10.5) 

GK,  GV,  Kl,  K2,  K3 

(5F10.5) 

Al,  A2,  A3,  Bl,  B2,  B3 

(6F10.5) 

NSAT 

(15) 

(BSTOP  (I),  BVELMAX  (I),I-1,  NV) 

(8F10.5) 

TABLE  2  -  PROGRAM  NOTATION  (INPUT) 


ANGLE 

AREA 

Al,  A2,  A3 
BM 

Bl,  B2,  B3 
BSTOP 
BVELMAX 
DCLDBFS 

DISPTON 

DUC 

GK 

GM 

GV 


Increment  of  wave  energy  spreading  for  calculation 
of  short-crested  responses;  e.g.,  5,  10,  or 
IS  degrees 

Fin  area,  square  feet  or  square  meters 
Fin  servo  coefficients 

Transverse  metacentric  radius,  (distance  between 
center  of  buoyancy  and  metacenter),  feet 
or  meters. 

Fin  controller  compensation  coefficients 

Fin  limit  angle 

Fin  limit  angular  velocity 

Free  stream  lift  coefficient  curve  slope, 
per  degree 

Ship  displacement,  long  tons  or  metric  tons 
Roll  damping  or  roll  damping  coefficients 
Overall  fin  control  gain 

Transverse  metacentric  height,  feet  or  meters 
Speed  dependent  fin  control  gain 


HO,  HI,  H2,  H3,  H4  Fin  lift  curve  correction  coefficients 


IPRINT 


ISC 

ITERATE 

IUNITS 
Kl,  K2,  K3 


L 


Array  of  print  options;  IPRINT(l)  *  SPECTRA, 
heading  printing  of  long-crested  response 
spectra  and  components.  IPRINT(2)  ■  ITERATN, 
step  by  step  printing  of  iteration  over  roll 
damping  of  relevant  variables. 

Switch  for  short-crested  responses.  If  ISC  + 

0,  provide  ANGLE. 

Equals  0  for  roll  damping  value  independent  of 
roll  angle.  Equals  1  when  iteration  is  required. 

Switch  indicating  type  of  units  for  Input/Output: 
IUNITS  ■  0  for  English,  -  1  for  metric. 

Roll  angle,  velocity,  and  acceleration  charac¬ 
teristic  gain  factors.  Sensitivities  of  demanded 
fin  angle  to  roll  angle,  velocity,  and 
acceleration. 

Ship  length  between  particulars,  feet  or  meters 
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TABLE  2  (Continued) 


M 

NAME1,  NAME 2,  NAME 3 

NSAT 

NSTAB 

NV 

NWH 

Q 

R 

SHAPE 

SWH 

T 

TITLE2 

TPHI 

TO 

VK 


Number  of  fin  pairs 

Identification  of  run  of  program  by  person's 
name,  code,  and  telephone  extension 

Equals  1  for  inclusion  of  significant  saturation 
effects 

Number  of  fin/controller/servo,  etc.,  input 
st  ts;  e.g. ,  £  10 

Number  of  ship  speeds,  e.g.,  5 

Kumber  of  sea  conditions,  e.g.,  _<  10 

6l/(I+6l),  see  Equation  (5)  for  definition 
of  6l/i 

Fin  moment  arm,  feet  or  meters 

Alphameric  description  of  waterplane  shape; 
e.g.,  PARAB,  ELLIP,  or  RECTANG 

Significant  wave  height,  feet  or  meters 

Draft,  feet  or  meters 

Alphameric  identification  of  a  run  of  the 
program;  e.g.,  ship  name 

Roll  period,  seconds 

Modal  wave  period,  seconds 

Ship  speed,  knots 


water  roll  decay  data,  or  by  use  of  analytically  predicted  values,  prior 
to  execution  of  FINCON.  The  iteration  is  performed  until 


|l-(a|/y)|  <  0.0!  (4) 


where  a?  is  that  RMS  roll  rate  obtained  from  the  prediction  using  the  roll 

<f> 

damping  coefficient  n  associated  with  y  by  Equation  (3). 


*In  practice,  a  straight  line  frequently  provides  an  adequate  repre¬ 
sentation  of  the  experimentally  derived  roll  data. 


Card  12  defines  two  print  options.  If  IPRINT(l)  is  SPECTRA,  a  heading 
by-heading  printing  of  long-crested  roll  spectra  and  their  components 
(e.g.,  the  wave  frequencies  ta,  the  response  amplitude  operators  (RAO's), 
etc.)  will  be  printed.  If  IPRINT(2)  is  ITERATN,  step-by-step  printing  of 
the  iteration  over  roll  damping  of  the  relevant  variables  will  occur.  This 
option  is  useful  for  debugging  purposes. 

If  NSTAB  is  zero  on  card  13,  only  unstabilized  roll  predictions  will 
be  made  and  cards  14  to  17  can  be  eliminated.  If  NSTAB  is  greater  than 
zero,  cards  14  to  17  should  be  repeated  NSTAB  times. 

Card  18  specifies  whether  or  not  saturation  calculations  will  be  done. 
If  NSAT  is  zero  (blank  card)  no  other  cards  are  needed.  If  NSAT  is  greater 
than  zero,  NSAT  pairs  of  values  for  fin  limit  angle  (BSTOP)  and  fin  limit 
velocity  (BVELMAX)  are  needed.  If  BVELMAX  is  not  supplied  (field  left 
blank)  a  value  will  be  generated  by  FINCON  (e.g.,  BVELMAX-10*BSTOP/TPHI) . 

Table  3  gives  a  sample  listing  of  input  cards  for  an  example  ship, 
between  the  END  OF  RECORD  and  END  OF  FILE  cards.  It  should  be  noted  that 
the  input  values  are  given  in  units  of  feet  and  long  tons.  A  metric  con¬ 
version  option  may  be  invoked  to  allow  the  input  of  values  in  metric  units. 
The  method  for  activating  this  option  is  to  place  a  "1"  in  Column  15 
(IUNITS)  of  card  3. 

The  first  two  data  cards  shown  in  Table  3  contain  only  alphanumeric, 
or  descriptive  data.  Card  3,  indicates  that  two  sea  conditions  and  one 
speed  are  to  be  considered  and  that  English  units  are  assumed  for  Input/ 
Output.  Cards  4  and  5  define  the  sea  conditions  in  terms  of  significant 
wave  height  and  modal  wave  period.  Card  6  specifies  the  ship  speed  at  25 
knots.  Card  7  gives  the  ship  particulars  of  displacement,  length,  draft, 
transverse  metacentric  height,  transverse  metacentric  radius,  roll  period, 
and  Q.  Q  is  the  ratio  of  added  mass  to  total  mass  moment  of  inertia, 
6l/(I+6l),  and  as  shown  in  Reference  1,  can  be  estimated  from 

-  0.186  +  1.179  Cg  -  0.615  C g2  (without  bilge  keels) 
and  (5) 

7^  -  -  0.002  +  0.814  C  -  0.316  C  2  (with  bilge  keels) 


TABLE  3  -  TYPICAL  CONTROL  AND  DATA  CARD  SET 


CHHJZXX .CM77000 ,T199 ,P4 . 

CHARGE, CHHJ.HQHAA15023.RS. 

ATT  ACH ,Ot DPL , F I NCONRE V I S I 0N6 , I D* PUAA . 
MAP (ON ) 

SETCORE(INDEF.ADDR) 

ATTACH(NSRDC) 

LIBRARY(NSRDC) 

FTN,I=OLDPL. 

LGO. 

7/8/9  END  OF  RECORD 


1  1: 

H. JONES 

1568 

71210 

2: 

CGN-42 

SELECTED  FINS  (TWO  FIN  PAIRS, 

75  SQ.FT.)  R.  N. 

3: 

2 

1  0 

4: 

24.61 

18.04 

5: 

12.9 

12.3 

6: 

25.0 

7: 

12000.0 

560.0 

22.7 

4.3 

16.06  12.8 

8: 

PARAB 

9: 

1 

15.0 

10: 

1 

11: 

0.1693 

0.00570 

-.0002 

12: 

13: 

1 

14: 

2 

75.0 

33.45 

15: 

0.43 

0.349 

1.117 

-0.519 

16: 

1.0 

1.650 

1.0 

2.5 

1.0 

17: 

1.0 

0.160 

0.025 

1.0 

0.630  0.092 

18: 

1 

19: 

10.2 

7.98875 

6/7/8/9  END  OF  FILE 


(25  KTS) 


0.331 


where  Cg  is  the  block  coefficient  and  average  bilge  keels  are  assumed. 

Card  8  specifies  waterplane  shape  as  parabolic,  which  is  usual  for  fine- 
form  naval  ships.  Card  9  indicates  that  short-crested  seas  will  be  treated 
and  specifies  the  spreading  angle  to  be  15  degrees.  If  card  9  were  blank, 
the  program  would  assume  that  only  long-crested  calculations  would  be  done. 
Card  10  indicates  that  iteration  over  the  roll-roll  damping  curve  is  re¬ 
quired  to  account  for  nonlinear  roll  damping.  Card  11  gives  the  coeffi¬ 
cients,  DUC  or  d,  of  the  roll-roll  damping  curve  (e.g.,  see  Equation  (3)). 
Had  roll  damping  been  linear  in  this  sample  input  case,  card  10  would  have 
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been  blank  and  card  11  would  have  contained  the  single  value  for  roll 
damping  coefficient,  n,  for  25  knots.  Card  12  is  blank,  so  printing  of  any 
intermediate  steps  (e.g.,  the  long-crested  spectra  and  their  RAO's,  etc., 
or  the  iteration  steps)  is  not  done. 

Cards  1  to  12  provide  all  the  data  necessary  to  complete  a  FINCON  run 
to  calculate  unstabilized  roll  angles.  If  card  13  were  blank,  this  is 
exactly  the  way  the  program  would  execute.  However,  card  13  indicates  that 
one  set  of  stabilizing  conditions  is  to  be  considered.  Cards  14  to  17 
provide  the  fin  and  fin  controller  particulars.  Two  pairs  of  75  square 
feet  fins  with  a  fin  moment  arm  of  33.45  feet  are  indicated  on  card  14. 

The  fin  moment  arm  is  taken  about  the  longitudinal  axis,  through  the  center 
of  gravity,  and  measured  to  the  center  of  the  fin  for  all  fin  pairs.  Card 
15  gives  the  free  stream,  lift  coefficient  curve  slope  and  the  lift  cor¬ 
rection  coefficients  which  compensate  for  fin- induced  sway  and  yaw  motions. 
The  values  of  card  15  are  estimated  by  the  procedures  outlined  in 

Reference  1.  Card  16  specifies  both  the  overall  gain  G  and  speed  dependent 

K. 

fin  control  gain  Gy,  as  well  as  the  roll  angle,  roll  rate,  and  roll  accel¬ 
eration  characteristic  gains  k^,  k2,  and  k^,  such  that  the  fin  controller 
equation  is 

B  =  Gk  •  G^kj*  +  k2i  +  k3V)  (6) 

where  <f>  is  the  roll  angle  and  6  is  fin  angle.  Card  17  specifies  the  fin 
servo  coefficients  and  the  fin-controller  compensation  coefficients.  The 
values  given  in  Table  3  for  card  17  variables  were  selected  from  Reference 
4  and  are  also  given  in  Reference  1. 

Card  18  is  not  blank,  indicating  that  any  significant  saturation 
effects  will  be  included  in  the  calculations,  and  so  an  additional  card  is 
required  to  provide  the  limiting  angle  and  limiting  speed*  of  the  fins. 
Insignificant  saturation  effects  are  automatically  ignored  by  the  program. 


*If  the  angle  is  provided  and  the  speed  is  left  blank,  the  program 
will  compute  an  appropriate  value. 
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This  is  to  avoid  Che  addicional  cosC  which  would  ocherwise  be  incurred  while 
having  no  significanC  effecC  on  Che  resulcs.  Had  card  18  been  blank,  indi- 
cacing  chac  no  sacuracion  effects  were  Co  be  considered,  Chen  no  ocher  daCa 
cards  would  have  been  needed. 

PROGRAM  OUTPUT 

Table  4  presents  Che  program  ouCpuC  for  Che  sample  inpuC  given  in 
Table  3.  The  firsc  page  ouCpuCS  Che  inpuC  idencifying  CiCles;  Che  second 
page  ouCpuCs  Che  operacing  condicions  and  ship/fin/fin  conCroller  particu- 
lars;  and  Che  Chird  page  provides  Che  results.  The  first  listings  on  the 
third  page  gives  the  resulting  unstabilized  RMS  roll  angles  and  correspond¬ 
ing  damping  coefficient  values  for  ship  headings  from  0  to  180  degrees 
(following  to  head  seas)  in  15-degree  increments  for  short-crested  seas. 

The  next  row  of  tables  contains  the  corresponding  RMS  stabilized  roll 
angles  and  damping  coefficients,  as  well  as  the  resulting  RMS  fin  angles 
and  velocities.  Those  values  of  RMS  roll  for  which  saturation  effects 
would  produce  less  than  a  2  percent  change  are  indicated  by  an  asterisk. 

For  such  headings,  only  the  unsaturated  values  are  calculated. 

Had  there  been  a  second  stabilized  condition  specified  on  input  card 
13,  another  row  of  tables  would  follow  the  one  for  case  1.  Results  for 
additional  speeds  and  sea  conditions  would  be  printed  in  a  similar  fashion 
on  subsequent  pages.  The  fourth  page  indicates  that  the  program  completed 
execution  satisfactorily;  e.g.,  no  system  (loader,  input/output,  etc.) 
errors  were  encountered. 

Table  5  presents  a  typical  output  when  IPRINT(l)  is  SPECTRA  on  card 
12.  One  such  page  would  appear  for  each  ship  heading.  The  columns  provide 
wave  frequency  W,  wave-encounter  frequency  WE,  wavelength  LAM,  wavelength- 
to-shiplength  LAM/L,  wave  number  K,  nondimens ional  transfer  function  TR, 
nondlmensional  response  amplitude  operator,  RAO,  wave-slope  spectrum  W  SL  S, 
and  roll  (unstabilized)  response  spectrum  SUR.  Also  given  are  the  dimen¬ 
sional  response  amplitude  operator  RAO  DIM,  the  wave-height  spectrum  W  HT  S, 
and  the  resulting  roll  response  spectrum  SUR  DIM,  which  should  be  equiva¬ 
lent  to  SUR.  Due  to  the  fact  that  this  sample  is  for  the  case  of  nonlinear 


TYPICAL  PROGRAM  OUTPUT,  ITERATION  OVER  ROLL  DAMPING 


*  *  ROLL  MOTION  PREDICTION  PROGRAM  •  *  • 


H. JONES 


1568 


71210 


!»19  O.ouo  0.000  1.000  I.fcSO  1.000  2.500  1.000  1.000  .100  .02$  1.000  .030  .092 


TABLE  4  (Continued) 


•••END*** 
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TABLE  5  -  TYPICAL  PROGRAM  OUTPUT  FOR  IPRINT(l)  =  SPECTRA  OPTION 


HU)a<*ao«nnw9«4NiB<4N«4iANKMKV4«na«ibiK#Mn(y<4f4«4Ci«o«« 


•4  IfN  «  *  <g  *4 
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Hi0O«49KlO  #r)NN«4*4 


»»«9<4|MBNM<4lB9>««4MB^nnNN<409» 

>»»CU'^/«W^OKl/S^#WMWWWK)WnMN 

HN#M^nM«4«<<»oeaee(»oaoo04« 


4O9>K>I0#NSI 


ienN#«t9N«4M«9M09in«4)My0IB^rM»N»O#9g»f4a«MB4f 

ianniBf*i0NiB«4N<HUSNK^^K«4tfna«M(t#nnMv««4Wo«»e( 

)eoo»)99iiii*<K«eNSiM»MM«4f4«4a««tas>0aao«ot( 


wu«aa0«0*09M'U|OoKi0i0ifitfiueeK»»tf> 

V)OOOOaH»OONKN^»i04nNM9«KN«4 


9  *  *4  ^  ISI  «  «  I 
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«  Ul 
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o  m 
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roll  damping,  the  spectral  data  and  the  RMS  roll  of  Table  5  are  not.  In 
general,  considered  especially  meaningful.  Had  the  damping  been  linear,  as 
Is  usually  the  case  for  the  ship  without  bilge  keels,  the  spectra  would  be 
representative  of  the  long-crested  seas  case  and  the  RMS  roll  could  corre¬ 
spond  to  a  value  given  on  the  "page-three-type-output"  of  Table  4. 

Table  6  presents  a  typical  output  when  IPRINT(2)  Is  XTERATN  on 
card  12.  In  brief,  the  intermediate  roll  rate  and  damping  coefficient 
values  are  printed  for  each  heading,  speed,  and  sea  condition.  Appendix  C 
further  describes  this  output  example. 

TABLE  6  -  TYPICAL  PROGRAM  OUTPUT  FOR  IPRINT(2)  =  ITERATN  OPTION 


IV 

IMU 

NTRY 

PHIN 

YP+1 

YP 

YP-1 

GP 

GP-1 

1 

7 

1 

.094 

.475 

0.000 

.239 

4.628 

.476 

1 

7 

1 

.094 

.475 

2.314 

0.000 

4.628 

4.628 

1 

7 

2 

.140 

.475 

2.314 

0.000 

3.847 

4.628 

1 

7 

2 

.140 

3.461 

3.461 

2.314 

3.847 

3.847 

1 

7 

3 

.156 

3.461 

3.461 

2.314 

3.674 

3.847 

1 

7 

3 

.156 

3.645 

3.645 

3.461 

3.674 

3.674 

1 

7 

4 

.159 

3.645 

3.645 

3.461 

3.652 

3.674 

1 

7 

1 

.094 

3.645 

0.000 

3.461 

.507 

3.674 

1 

7 

1 

.094 

3.645 

.254 

0.000 

.507 

.507 

1 

7 

2 

.100 

3.645 

.254 

0.000 

.506 

.507 

1 

7 

2 

.100 

.504 

.504 

.254 

.506 

.506 

] 

7 

3 

.  105 

.504 

.504 

.254 

.504 

.506 

PROGRAM  EXECUTION 

A  typical  deck  control  card  set  up  is  given  in  Table  3.  Simply 
speaking,  the  object  program  is  attached  and  executed.  The  object  program 
is  stored  permanently  on  a  private  disk  pack  and  can  be  recovered  for 
storage  on  the  main  (public)  disk  and  for  user  execution  by  running  the 
control  card  deck  of  Table  7.  The  source  deck  is  also  stored  on  the  pri¬ 
vate  disk  pack  in  an  UPDATE  file  such  that  program  modifications  can  be 
easily  made,  if  necessary.  The  program  listing  of  Appendix  B  was  printed 
from  this  UPDATE  file. 
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TABLE  7  -  CONTROL  CARD  SET  TO  RETRIEVE  OBJECT  PROGRAM 


COU.  123456789112345678921234567893123456789412345678951234567896123456789 

CHHJPAK,CM77777,T100,RP1,P3. 

CHARGE, CHHJ,XXXXXXXXXX,CC,R. 

PAUSE.  JOB  REQUIRES  DISK  PACK  DV4850. 

MOUNT ,VSN=DV4850 ,SN=HJPKL4 . 

REQUEST, TWO, *PF. 

ATTACH .ONE ,F INCONOBJECTNOV , I D=CHHJ ,CY=1 ,MR=1 , SN=HJPKL4 . 
COPYBF.ONE.TWO.l. 

CATALOG, TWO, F INCONOBJECTNOV, ID=PUAA,AC=XXXXXXXXXX,CY=1,MR=1. 

6/7/8/9  END  OF  FILE 


The  run  time  of  the  program,  indicated  by  TXXX  on  the  job  card  of 
Table  3,  varies,  of  course,  with  the  amount  of  calculation  required. 
Roughly,  for  nonlinear  roll  damping  and  15-degree  short-crested  spreading, 
a  stabilized  roll  calculation  (without  saturation  effects,  and  for  a  single 
sea  condition  and  speed)  takes  about  50  seconds  of  execution  time  and  15 
seconds  compilation  time.  For  unstabilized  calculations,  the  execution 
time  is  somewhat  less  than  half  of  the  time  for  the  stabilized  case.  The 
time  increases  proportionately  as  the  spreading  angle  of  the  short-crested 
seas  is  decreased.  For  a  5-degree  spreading  angle,  the  time  is  almost 
three  times  145  seconds)  that  of  the  15-degree  case  (reflecting  the 
fact  that  there  are  about  three  times  as  many  calculations  that  need  to  be 
performed).  From  the  runs  made  to  date,  it  is  not  evident  that  decreasing 
the  spreading  angle  from  15  degrees  increases  the  accuracy  of  predicted 
roll  at  a  given  ship  heading  Ky  a  noticeable  amount.  However,  a  finer  mesh 
of  spreading  angles  does,  in  some  cases,  permit  a  more  refined  localization 
of  the  worst  heading  angle.  Thus,  the  required  execution  time  is  a 
multiple  of  50  seconds  depending  on  the  number  of  speeds,  sea  conditions, 
and  the  value  of  the  spreading  angle  (in  proportion  to  15  degrees). 
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The  required  memory,  as  indicated  by  CMXXXXXX  on  the  job  is  mil 
octal  words  (see  Table  3).  The  job  priority,  indicated  by  PX  on  the  job 
card,  is  then  determined  by  the  amount  of  system  time  required.  Based  on 
current  computer  center  figures  for  the  CDC  6700  and  average  costs  over 
several  program  runs,  the  guidelines  in  Table  8  are  offered. 

TABLE  8  -  RUN  TIME  AND  COST  GUIDELINES 


T  Highest  Priority  (Turnaround) 

<  200  4  (prime  shift,  1  hour  max  after 

completion) 

<3600  3  (prime  shift,  as  soon  as  possible) 

Unlimited  2  (nonprime  shift,  when  possible, 

overnight) 


$/ System  Seconds 
0.090 


0.074 

0.060 


PROGRAM  VERIFICATION 

Predicted  values  of  unstabilized  ship  roll  motion  using  the  one- 
degree-of -freedom  roll  motion  procedure,  have  been  compared  to  model  and 
full-scale  experiment  results  in  both  References  1  and  2.  Ongoing  work  at 
DTNSRDC  by  Meyers  has  found  the  results  of  the  single-degree  equation 
very  similar  to  those  of  the  coupled,  three-degree  equations  for  roll-sway- 
yaw  for  the  worst  heading  roll  motion,  although  some  underprediction  in  bow 
seas  and  some  overprediction  in  following  to  quartering  seas  have  been 
noted.  Additionally,  Reference  2,  as  well  as  work  by  Lloyd  and  other 
Admiralty  Marine  Technology  Establishment  (AMTE)  personnel,  has  substan¬ 
tiated,  at  least  in  part,  the  stabilized  roll  and  fin  predictions.  It  is 
generally  recognized  that  the  predictions  of  FINCON  (e.g.,  for  the  worst 
heading)  are  appropriate  for  use  in  design  problems. 

The  coding  of  FINCON  has  been  verified  by  making  comparisons  with 
results  of  the  older  unpublished  FINS  program,  as  well  as  with  published 
results  of  programs  currently  used  by  Lloyd  and  others  at  AMTE.  The  com¬ 
parisons  substantiate  the  correctness  of  the  coding  of  FINCON  in  general, 
though  some  differences  do  occur  between  the  results  of  FINCON  and  the 


AMTE  program.  For  example,  at  some  speeds,  the  FINCON  unstabilized  roll 
angles  were  higher  than  the  corresponding  AMTE  values  for  the  same  narrow- 
beam  LEANDER-class  frigate  evaluated  for  the  case  when  roll  damping  is  in¬ 
dependent  of  roll  angle.  The  differences  may  be  due  to  differences  in  the 
input,  slight  differences  in  the  motion  equations  and  algorithms  programmed, 
or  differences  in  the  short-crested  seas  algorithm.  One  known  difference 
is  that  FINCON  accepts  as  input  the  true  value  for  BM,  the  center  of 
buoyancy;  whereas  the  AMTE  program  computes  a  value  based  on  the  shape  of 
the  waterplane.  Another  difference,  though  not  relevant  to  the  comparisons 
made  for  the  LEANDER,  is  that  the  AMTE  program  has  no  provision  to  handle 
the  case  when  roll  damping  is  dependent  on  roll  angle. 

FUTURE  WORK 

FINCON  is  the  fJrst  of  a  series  of  new  tools  being  developed  to 
enhance  the  U.S.  Navy's  roll/fin  design  capability.  As  such,  FINCON  is 
the  basis  for  all  such  current  and  near-future  investigations.  Specific 
guidelines  for  optimum  use  of  the  program  in  the  form  of  a  rather  complete 
design  exercise  can  be  found  in  Reference  1.  Procedures  for  evaluating 
bilge-keel  and  fin  sizing  and  stabilizer  control  optimization  are  detailed 
there.  These  procedures  indicate  how  the  use  of  the  FINCON  program  in  roll/ 
fin/controller  design  practice  can  be  extremely  instructive. 

Another  very  important  area  currently  being  investigated  by  Cox  is  the 
use  of  a  coupled,  three-degree  system  of  equations  for  roll,  sway,  and  yaw 
motion  prediction.  A  more  general  and  refined  program  is  being  developed 
in  conjunction  with  that  work.  The  program  under  development  will  also 
be  of  practical  use  at  an  early  stage  of  ship  design,  requiring  only  very 
simple  input  requirements.  A  complete  report  and  user's  manual  for  the 
improved  three-degree-of-freedom  simulation  system  will  soon  be  published. 

It  will  include  details  of  the  approach  which  is  used  in  the  current  one- 
degree-of -freedom  program  to  recognize  fin-saturation  effects. 

CONCLUDING  REMARKS 

This  report  provides  a  user's  guide  to  FINCON,  a  roll,  fin,  fin  con¬ 
troller  prediction  computer  program.  No  attempt  to  describe  design 
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practices  or  the  required  engineering  decisions  necessary  to  using  this 
tool  has  been  made;  Reference  1  provides  a  comprehensive  discussion  of  such 
materials.  Sample  inputs  and  outputs,  as  well  as  a  description  of  the 
program  organization  and  procedures  have  been  given.  It  is  envisioned 
that  the  engineer,  with  a  working  knowledge  of  Reference  1,  will  run  the 
program  essentially  as  a  "black  box" — he/she  is  not  expected  to  need  to 
contend  with  the  actual  FORTRAN  or  source  deck;  and,  thus,  only  a  very 
rudimentary  knowledge  of  programing  or  computers  is  required.  Instead, 
he/she  will  be  required  to  actively  participate  in  the  engineering  tradeoff 
decisions  necessary  in  design  work,  and,  as  such  will  probably  run  FINCON 
several  times  in  any  given  investigation. 


APPENDIX  A 

PROGRAM  PROCEDURE  AND  FLOW 


The  equations  solved  in  FINCON  are  listed  In  Table  9  and  are  taken 
almost  exclusively  from  Reference  1.  For  purposes  of  Illustration,  a 
short-crested  spreading  angle  of  10  degrees  Is  assumed.  Table  10  provides 
a  description  of  the  nomenclature  used  in  Table  9.  The  corresponding 
FORTRAN  notation  (e.g.,  see  the  listing  of  Appendix  B)  follows  as  closely 
as  possible  that  of  Table  10.* 

Table  9  is,  in  a  sense,  broken  into  three  algorithmic  steps.  Four 
basic  predictions  are  identified:  stabilized  roll,  roll  rate,  fin  angle, 
and  fin  velocity.  Each  of  the  four  is  identified  by  a  more  or  less  reverse 
building-block  procedure.  For  example,  the  final  step  of  the  first  al¬ 
gorithm  is  labeled  1-1,  the  step  preceding  1-1  is  1-2,  the  step  preceding 
1-2  is  1-3,  etc.  Similarly,  steps  2-1,  2-2,  etc.,  and  3-1,  3-2,  etc.,  are 
developed.  It  is  felt  that  this  reverse  building-block  approach  to  listing 
the  steps  makes  it  easier  to  see  the  final  results  and  is  also  representa¬ 
tive  of  the  procedure  followed  in  organizing  the  equations  of  Reference  1 
for  programming  purposes.  Obviously  some  of  the  steps  developed  for  the 
first  algorithm  are  needed  by  the  other  two  algorithms  (e.g.,  step  1-4-1); 
however,  it  was  not  considered  necessary  to  repeat  these  for  each  of  the 
other  two.  Instead,  for  clarity,  one  can  assume  that  the  results  of  each 
step  of  the  first  algorithm  are  available  to  the  remaining  algorithms. 

Figure  3  presents  a  diagram  of  the  flow  sequence  of  FINCON.  The 
figure  identifies  the  important  loops  over  sea  conditions  and  ship  speed 
for  both  unstabilized  and  stabilized  predictions.  The  diagram  was  con¬ 
structed  with  the  intent  of  providing  a  quick  overview  of  the  entire 
program  flow  so  that  major  computational  segments  are  easily  identified. 


*0ne  exception  to  this  is  that  c 
FORTRAN.  u 


in  Table  9  becomes  CA  in  the 


TABLE  9  (Continued) 


TABLE  9  (Continued) 


Ge  *  71  "VW1  -  P«' 

L 


Sin  kL*  1  * 

C  =  — — ;  1^*  *  j  kL  cos  v» 

L 

L  -  L,  1/2  L,  or  /7/4  L  for  rectangular,  parabolic,  or 
elliptical  waterplane,  respectively; 


1-4-1 


i  \  2  /  v  2 

M  /CU 


1+2 


cos  5  +  ^  Sin  k)  +  (^f- 

S  /V  8  S 


1-4-2  c  *  (a2+b2)1/2;  b  *  2n  A 
-  s  s  s  s 


(pv2) 


•  -  -  'fin  1 


—  —  -  c  <j>  ^sm  AGM 

8  8 


(?) 


<®a>o 


,  2  2.1/2  4> 
Cs(aR+al> 


1-4-4 


sin  £ 


kI(aRbR~aIbI)*kR(aRbl'>'aIbR) . 

2  2.  „2,,  2.  ,1/2  ’ 

[(kR+ki)(aR+ai)(bR+bi)] 


cos  £ 


tR(»RbR-,IbI>',~tI(,Rbl't-*Ib8) 
.J,,!,  2,  2.  ,,2. .2.  .1/2 
l(kR*kI)(,R+*I)<bR+1’l)) 
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TABLE  9  (Continued) 


TABLE  10  -  NOTATION  USED  IN  EQUATIONS  OF  TABLE  9 


A 

Fin  area 

a 

1-  A2 

aI 

V2 

* a2 * 

Fin  servo  coefficients 

BM 

Distance  between  center  of  buoyancy  and  metacenter 

bi 

“Eb2 

bR 

»1  -  UEb3 

Fin  controller  compensation  coefficients 

(aW  )1/2 

u,s 

c 

u,s 

(dCL/d3)E 

Effective  lift  curve  slope 

<dCL/dp)FS 

Free  stream  lift  coefficient  curve  slope 

GM 

Transverse  metacentric  height 

Gk 

Fin  controller  overall  gain  control 

Gv 

Fin  controller  velocity  dependent  gain  control 

g 

Acceleration  due  to  gravity 

W 

Ratio  of  effective  to  free  stream  fin  lift  curve  slopes 

ho,h^,h2»h2»h^ 

W  "  ho  +  +  h2a)E  +  h3WE  +  h4WE 

Jo,Jl 

Bessel  functions  of  the  first  kind 

k 

2 

Wave  number,  to  /g 

kI 

“Ek2 

kl  '  “Ek3 

30 


TABLE  10  (Continued) 


k^  1 1^2  > 

L 

H 

n 

u 

q 

R 


B(Ba) 

estp 

^max 

A 


?a 

(Kh/3 

A 


X 


Roll  angle,  velocity,  and  acceleration  characteristic 
gain  values 

Ship  length  between  perpendiculars 
Number  of  fin  pairs 

Roll  decay  coefficient  (ship  without  fins) 

6l/(I+6l),  see  Equation  (5)  for  definition  of  61/1 

Fin  moment  arm  (about  a  longitudinal  line  through  the 
ship  center  of  gravity) 

Saturation  multiplier 

Wave-slope  spectral  coordinate 

Unstabilized  roll  response 

Mean  ship  draft 

Modal  (peak)  period  of  wave-height  spectrum 
Unstabilized  roll  transfer  function 

Ship  speed 

Fin  angle  (amplitude) 

Fin  limit  angle 

Fin  limit  velocity 
Ship  displacement  weight 
Wave  amplitude 

Significant  wave  height 
Tuning  factor, 

Wavelength 

Ship  heading,  with  respect  to  the  ship,  of  predominant 
wave  direction 

31 


•- - 'full  IN  ifiTiftfc 


TABLE  10  (Continued) 


Wave  direction,  with  respect  to  the  ship,  apart  from 
the  predominant  direction 

Density  of  seawater 

Stabilized  root  mean  square  roll  angle 
Root  mean  square  fin  angle 
Root  mean  square  fin  velocity 
Stabilized  roll  angle  amplitude 
Unstabilized  roll  angle  amplitude 
Roll  reduction  factor 
Wave  frequency 
Frequency  of  wave  encounter 

Ship  natural  roll  frequency 

Frequency  above  which  roll  response  is  negligible, 

2  radians/second 


Long-crested 

Short-crested 

Stabilized 


Unstabilized 


(Sea  Conditions) 


FINCON  PROGRAM  FLOW 
I/O 

Initialize  , 

Calculate  cos^  values 
DO  700  (Sea  Conditions) 

Calculate  2-Param.  Wave  Slope  Spectrum 
r~  DO  600  (Ship  Speed) 

r-  DO  500  (Predominant  Wave  Direction) 

Initialize 

80  Next  unstablllzed  roll  damping  value 
_  DO  400  (Spreading  Angle)  — - 

—  Initialize 

—  r-  00  200  (Wave  Frequency) 

§  ~  Initialize 

«  u  Calculate  values  dependent  on  wave  frequency 

f  ||  Calculate  unstablllzed  roll  transfer  function 

?  o  Calculate  unstablllzed  roll  spectrum  (long  crested  only) 

3  200  End-of-Loop 

a  Calculate  unstablllzed  RMS  roll 

If  no  iteration,  store  long  crested  roll  values 

*—  400  End-of-Loop 

Sum  short  crested  components 
Iteration  over  roll  damping  required? 

No  -  go  to  step  250 

Yes  -  If  number  of  attempts  >  10,  write  message  and  go  to  step  900 

-  If  closed  go  to  step  250 

-  otherwise  perform  Iteration  procedure 

"g  and  go  to  step  80  • —  - 

—  If  requested,  output  long  crested  spectrum  and  Its  components 

u  250  If  stabilized  roll  calculations  not  required  go  to  step  500 

~  £  —  00  280  (Stabilizing  Cases) 

<>  5  80*  Next  stabilized  roll  damping  value 

o.  «  »—  DO  400*  (spreading  Angle) 

5  »-  00  200*  (Wave  Frequency) 

2  Z  _  Initialize 

vi  c  Calculate  effective  lift  curve  slope 

_c  u  -  tic  Calculate  fin  to  roll  angle  amplitude 

J—  cm3  Calculate  roll  reduction  factor 

>  ^  2  J  Calculate  stabilized  roll,  fin  motion,  and  fin 

>-  ecu.  velocity  spectra 

IT  -o  200*  End-of-Loop 

•  3  Calculate  stabilized  roll,  roll  rate,  fin  angle,  and  fin 

—  a  velocity 

If  no  iteration,  store  long  crested  roll  values 
—  Sum  Short  crested  components.  If  requested 

o  •—  400*  End-of-Loop 

“  Iteration  over  roll  damping  required? 

"S  No-go  to  250* 

■-  Yes  -  If  number  of  attempts  >  10,  write  message  and  go  to 

-  step  900 

m  "If  closed  go  to  step  250* 

v>  -  Otherwise  perform  iteration  procedure  and 

go  to  step  80* 

250*  No  (further)  damping  iteration  required 
L-  280  End-of-Loop 
L-  500  End-of-Loop 
Output  results 
—  600  End-of-Loop 
700  End-of-Loop 
900  End  of  run  message 

Note:  RMS  -  root  mean  square. 

Figure  3  -  Diagram  of  Program  Flow  Sequence 


(Unstabi lized  Roll) 


APPENDIX  B 
PROGRAM  LISTING 

A  complete  listing  of  FINCON  Is  given  on  the  subsequent  pages.  The 
routines  are  listed  in  the  order  FINCON,  ITREQ,  BRWSSP,  ALGRNG ,  FINSTAB,  and 
ERRF.  BESSJ  is  a  system  routine  and  hence  not  given  here.  The  listing, 
made  from  an  UPDATE  file  stored  on  a  private  disk  pack,  contains  a  unique 
identification  of  each  line  of  coding  on  the  far  right-hand  side  of  the 
pages.  This  identification  first  identifies  the  routine  by  a  name  (e.g., 
ROLL,  TREQ,  BUSS,  ALGR,  and  FST)  and  then  by  a  line  number. 

The  FORTRAN  is  embedded  with  comment  cards  for  identification  of  the 
steps  of  the  prediction  procedure.  Modifications  can  easily  be  made  to 
this  source  coding  by  inserting  or  deleting  statements  anywhere  in  the 
routines  via  an  "UPDATE"  run  of  the  program. 


PROGftA*  FXNCON 


74/74  OPT.Q  POUND* •/  TRACE 


FT*  4. *44*0 


11/26/7)  11.25.10 


1 


5 


10 


15 


70 


25 


30 


35 


40 


45 


50 


55 


PROGRAM  FIWCON  <INPUT*512.0UTPUT«512,TAPE5*IN»UT.TAPC6»OUTPUT» 

2  TAPED 

9 

•cnc  6700  -  -  -  MATCH,  1976  -  -  -  OTNSROC  -  -  -  COOE  1560,  SUSAN  BAlES 

•  J.  9.  TUCKER  •  -  -  AUGUST,  1979  -  •  CHI  ASSOCIATES,  INC.,  90SLYN,  VA, 
•FORTRAN  PROGRAM  TO  PREDICT  ROLL  NOTION  . • ,. UNSTABILIZE 0  NOTION  PREOICTE 
•USING  THE  THEORY  OF  J.  £.  CONOLLY  MODIFIED  TO  ALLON  FOR  NONLINEAR 

♦damping.  stabilized  notion  predicted  using  cox  and  lloyo. 

* 

C  THIS  VERSION  OF  FINCON  (TPSC76I  NRITES  TO  PF  THE  SMORTCRESTEO  STAS. 

C  OR  UNSTAS.  ROLL  RESPONSES  (RMSI.  THIS  FILE  IS  DESIGNED  TO  BE  POST 
C  PROCESSED  for  use  WITH  THE  polar  plot  program. 

C 

INTEGER  SHAPE .PARAS. ELLIP.RECTANG, SPECTRA 

REAL  KL,KL9RIN£,L,  NUR.MUR  ,  HOUR,  LPRI  ME  ,LAM,tAMCNL,Kl»K2,KJ,M 
CCMMON/ITRTN/STAT.PRCM.NTRY , KHU.KNU.KV  « GP.GPMl , YP91 ,YP, YPP1 « 

2  PHIN.IPSINTISt, ITERATE, NPHI 

COMMON  /STAB/  NSTAS.MdOl  ,  ARE  AIIOI.RTIOI,  OCLOBFS  (101,  HO  1101. 

2  Hllldt  ,H2C10t  ,M3(1CI  ,H4(10l,GK(10t  ,GV(1SI,K1(1S>  ,K2  1 101  ,KSU0  I  , 

2  A1(1O(,A2(1D>,A3(1OI,B1(10I,B2(1O(,B3(1O(,OANFU(13(» DAMPS (10,131. 
2  DUC  15,61  .SSIGLC  (10.1  31  .NNO.ILC. ISC, NE  11.0,351  , 

2  SSIGVLCdQ«13),SSIGVSC(10,13J, 

2  SUT(4O,35I.X(40I.NN,COSl(3«l,CONl,SSIGSCdO,13>, 

2  VFS (5>,DIS»LB.GH  ,  BMOTLC 1 10, 131 . 

2  BMOTSCdO,  13I,BVELSC110,13),  SVELLC  (10, 131, BST  OP  (51  •  SVEL  MAX  (5  I 
2  ,HSAT,ITEST(10.131 

DIMENSION  TITL£2(8I,  EXT  (4  0  l,P(351.S«»OI,SNHlt01,TO(101*VK(SI, 

2  HUR  (131  ,HUO(l  31  , NUR(  3  51 ,  SINNU  (351 ,  COSNU  (35)  ,  SIGSOSG  ( 13»  ,  SIGSC  (131 
2, SIGLC1 13I.TR (40 (.HOUR (351 .BESSEL ( 1001 , SURVI 40,35) 

2  «SGVSQSC  (1 3  )  ,  SIGVSC  (131  ,VH0UR(35I,  SIGVLCd  31  .SMMET  ( 101 
DATA  SPFCTRA , I TERATN/7HSPECTR A . 7HI TERATN/ 

DATA  MUO/0,15 ,30,45.60,75,90, 105,120,135.150,165,150/ 

OAT  A  PI ,RH3,  GPAVITT/3 .1415936 ,1.99, 32.1725/ 

DATA  PARA'S, ELLIP.RECTANG/  10NPARAB  .10HELLIP  . 

210HRECTANG  / 

DATA  NM/40/  ,  NMU/13/  .EPS/. 0001/,  NOEL/. 05/ 

• 

•INPUT  AND  OUTPUT  THE  SEA  AND  SHIP  CONDITIONS  FOR  UNSTA9ILIZEO  ROU 
•CALCULATIONS. 

• 

READ  <5.1000)  NANE1.NAHE2. NAMES 
WRITE  <6.2000)  NAHE1.MAME2. NAMES 

RFAO  <5,1000)  TITLE2 
WRITE  <6.2001)  TITLE? 

READ  <5.1001)  NHH.NV.IUNITS 

NRtTE  (1)  NWH.Ntf 

READ  <5.1002)  ISMMI1 ) .1  *  1 .MW Hi 

* 

IF  < I  UNI TS  .EO.  0)  GO  TO  5 
WRITE  <6.30021  <SNH<II.I*1,NWH) 

GO  TO  6 

9 

5  WRITE  <6.2002)  <$WM<I) ,I*1.NWH) 

6  continue 

RCAO  <5.10021  <TOUI,I*l.WMI 
WRITE  <t.?003)  <T0<I) .1*1 .NWH) 

READ  (5.10021  (VK(I).I*1.NV) 


ROLL 

2 

ROLL 

3 

ROLL 

4 

ROLL 

5 

ROLL 

6 

ROLL 

r 

ROLL 

o 

ROLL 

9 

ROLL 

IQ 

ROLL 

11 

ROLL 

12 

ROLL 

13 

ROLL 

14 

ROLL 

15 

ROLL 

IS 

ROLL 

17 

ROLL 

18 

ROLL 

15 

ROLL 

20 

ROLL 

21 

ROLL 

22 

ROLL 

23 

ROLL 

24 

ROLL 

25 

ROLL 

26 

ROLL 

27 

ROLL 

26 

ROLL 

29 

ROLL 

30 

ROLL 

31 

ROLL 

32 

ROLL 

33 

ROLL 

34 

ROLL 

35 

ROLL 

36 

ROLL 

37 

ROLL 

30 

ROLL 

39 

ROLL 

40 

ROLL 

41 

ROLL 

42 

ROLL 

43 

ROLL 

44 

ROLL 

45 

ROLL 

46 

ROLL 

47 

ROLL 

40 

ROLL 

49 

ROLL 

50 

ROLL 

51 

ROLL 

52 

ROLL 

53 

ROLL 

54 

ROLL 

55 

ROLL 

56 

ROLL 

57 

ROLL 

50 
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PROGRAM  FINCON 


74/74  OPT«0  ROUND. •/  TRACE 


F TN  4.6.660 


I1/M/7T  11. EG.  SO 


60 


61 


70 


75 


•0 


05 


40 


55 


100 


105 


110 


NRITE  (6*70001  ( VK  ( 1 1  ,I>1 ,  NVI  ROLL 

REAO  15.10071  OISPTON.L.T ,GH,BN,TPHI .0  ROLL 

•  ROLL 

IF  ( IUNITS  *EO.  01  GO  TO  7  ROLL 

MRITE  16,30061  OISPTOM.L ,  T.GM.BM.TPHI.O  ROLL 

GO  TO  4  ROLL 

•  ROLL 

7  MRITE  <6,70041  OISPTOM.L  ,  T  ,GN  .Bit,  TPHI  ,0  ROLL 

0  CONTINUE  ROLL 

RFAO  (5,10001  SHAPE  ROLL 

MRITE  (6,70051  SHAPE  ROLL 

REAO  (5,10061  ISC,  ANGLE  ROLL 

iFdsc.EQ.ot  go  to  io  Roll 

NRITE  (6,70061  ANGLE  ROLL 

10  REAO  (5.10061  ITERATE  ROLL 

STAT  .1.0  ROLL 

POOH  •  0.01  ROLL 

IF  (ITERATE  ■ GT .  01  MRITE  (6,7017)  ROLL 

IF  (ITERATE  .E9.0)  REAO  (5,10071  ( DUC ( IV.lt , IV« 1 ,NVt  ROLL 

IF  (ITERATE  •  EQ ,0 1  MRITE  (6,70161  (OUC (I», It , IV»1,NVI  ROLL 

IF  (ITERATE  .  EQ.O I  GO  TO  70  ROLL 

MRITE  <6.70071  ROLL 

00  15  1VM.NV  ROLL 

REA0<5, 10071  (DUCdV.lt  ,I»t,6t  ROLL 

15  MRITE<6. 71041  VKtlVt  , (DUCtIV.lt, 1*1, 6t  ROLL 

70  REAO  (5.100  01  IPRINT  ROLL 

1F<  < ISC .ME, 0 1 • ANO. ( I PRINT (It. EQ. SPECTRA 1 1  GO  TO  77  ROLL 

IF  ( IPRINTI 1 1  .EQ.  SPECTRA)  MRITE  (6,70101  ROLL 

IF  (IPRIMTtE)  ■  EQ •  HERAT N)  MRITE  (6,7014)  ROLL 

GO  TO  75  ROLL 

77  MRITE(6. 70741  ROLL 

GO  TO  400  ROLL 

•  ROLL 

•INPUT  ANO  OUTPUT  CONDITIONS  FOR  FIN  STAOILI7EO  ROLL  PREDICTIONS.  ROLL 

•  ROLL 

75  REAO  (5,10011  NSTIO  ROLL 

IF  (NSTAN  ,LT.  It  GO  TO  35  ROLL 

IF  (NSTAB  ,G£.  1  .ANO.  IUNITS  .EQ.  0)  MRITE  (6,70761  NSTAB  ROLL 

IF  (NSTAB  .JE.  1  .AND.  IUNITS  .EQ.  It  MRITE  (6,30761  NSTAB  ROLL 

OO  76  I«l, NSTAB  ROLL 

REAO  (5,10061  H(It.AREA(It,R(  It  ROLL 

REAO  (5.10071  OCLDBFS  (I)  ,  H0(I  I  ,H1  <  I  ),H7  (I  t.  H3(I I  ,  H6<  1 1  ROLL 

REAO  (5,10071  GK  ( 1 1  ,GV(I  I  ,  Kit  I)«K7dl  ,53  ( I)  ROLL 

REAO  <5,10171  Al(I),A7(t),A3(I),Bl(II,B7(II  .B3II)  ROLL 

70  MRITE  (6,70751  I ,  M(I )  .AREA  (It  ,  R 1 1 1 ,  DCLDBFSd  t  ,H0  (It  ,  Nidi  ,H7<  It  ,  ROLL 

7  H3dl.H6(II,GK(II,GV(It,Kl(I)  ,K7(I)  ,K3(II  ,Al(I),A7(II,A3(It,  ROLL 

7  Bit  It .87(11 ,83(11  ROLL 

REAO  (5,10011  NSAT  ROLL 

IF  (NSAT. £0.0 1  GO  TO  30  ROLL 

REAO  (5,10071  <B STOP (IV) .BVELNAXII Vt ,IV  •  l.NVt  ROLL 

DO  74  IV  >  l.NV  ROLL 

75  IF  (BVELNAK  IV)  .EO.O.Ot  BVELHAAlIVt  »  10. •BSTOPdVt /TPHI  ROLL 

MRITE  (6,70701  (IV.BSTOPdVt  .BVELNAX  (IVt  ,IV  *  l.NVt  ROLL 

30  CONTINUE  ROLL 

35  CONTINUE  ROLL 

•  ROLL 

•INITIALIZE.  ROLL 


54 

60 

61 

67 

63 
66 

65 

66 
67 
60 

64 

70 

71 
77 

73 
76 

75 

76 

77 
70 

74 
00 
01 
07 
03 
06 
OS 
06 
07 
00 
04 

40 

41 
57 

43 
54 

45 

46 

47 
40 

44 
100 
101 
107 

103 
106 

105 

106 
107 
100 

104 
111 
111 
117 
111 

114 

115 
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P“06*AM  FINCON 


24/7.  0PT*C  ROUND* •/  TRACE 


FTN  k. 6*1.60 


ll/H/Tf  tt.2f.ll 


115 

• 

ROLL 

in 

IF  1 1 UN ITS  . EQ.  01  GO  TO  30 

ROLL 

112 

• 

ROLL 

tt. 

IJH  =  3.20054 

ROLL 

tt« 

OO  36  I  «  l.NWH 

ROLL 

121 

120 

SMNcTIII  x  SWNC II 

ROLL 

121 

SMNIII  »  54*4  1  1 1  *UH 

ROLL 

122 

36 

CONTINUE 

ROLL 

121 

OISPTON  *  3ISPTON*0.9642 

ROLL 

12. 

L  *  l*UM 

ROLL 

12* 

125 

T  *  T *UN 

ROLL 

12* 

G*  =  G**UH 

ROLL 

127 

ON  *  ON*UN 

ROLL 

12* 

• 

ROLL 

12* 

IF  INSTA9  •  L T .  11  GO  TO  30 

ROLL 

til 

130 

00  37  I  *  1 . NSTAB 

ROLL 

lit 

ARE  Alii  *  ARE  A 1 X) •10*7639 

ROLL 

111 

Rfll  *  Rm*u« 

ROLL 

113 

37 

CONTINUE 

ROLL 

tl. 

• 

ROLL 

11* 

135 

35 

CONTINUE 

ROLL 

11* 

• 

ROLL 

112 

MPNI*?.*PI/TPHI 

ROLL 

13* 

0IS»L9  =  2240.  •  OISPTON 

ROLL 

13. 

00  -0  l¥rl,NV 

ROLL 

1.1 

140 

4  0 

VFSIIVI  X  1.6576*  VKIIVI 

ROLL 

1.1 

TE^Nt  *  2.*PI*GRAVITV/L 

ROLL 

1.2 

00  60  IH*1*NN 

ROLL 

1.1 

NttN  1  xFLOAT  UNI  *NOE\. 

ROLL 

1.. 

60 

PNTIINI  x  £*P  I-W(IM) ••2/GRA¥I TV  •  Tl 

ROLL 

1.1 

145 

• 

ROLL 

1.6 

•GENERATE  VALUES  r0R  COSINE  SQUARED  SPREADING  OF  NAVE  ENERGY. 

ROLL 

1*2 

• 

ROLL 

1.1 

IF  ( ISC  .£3.  01  GC  TO  75 

ROLL 

1.. 

HA  x  (lftn,'AN6LE>/2. 

ROLL 

1*1 

150 

CONI  *  l./NA 

ROLL 

l.t 

NNU  x  2*IFIR 1  HA  1  -  l 

ROLL 

1*2 

ILC  x  NNU/2  ♦  l 

ROLL 

1*1 

CON 2  x  PI/<2 . *HA> 

ROLL 

15. 

J  x  -IFIKCNAI 

ROLL 

1** 

155 

00  ."0  I*1*NNU 

ROLL 

1*6 

J  *  J  ♦  1 

ROLL 

152 

Pfll  x  J 

ROLL 

191 

70 

COSHII  *  COS  TP  III *CON?>  **2 

ROLL 

IS* 

75 

CONTINUE 

ROLL 

16* 

160 

* 

ROLL 

161 

•BEGIN 

LOOP  OVER  SEA  CONDITIONS. 

ROLL 

162 

• 

ROLL 

163 

00  7  0D  INH*1 « NNH 

ROLL 

16. 

• 

ROLL 

165 

165 

•COMPUTE  RRcTSCHNEIOER  2-PARAMETER  NAVE  SLOPE  spectrum. 

ROLL 

166 

• 

ROLL 

162 

CALL  RRHSS>  INN.SNH(INHI,TO(INHI.N,S) 

ROLL 

16* 

* 

ROLL 

16. 

•BEGIN 

LOOP  OVER  SNIP.  SPEED. 

ROLL 

121 

170 

9 

ROLL 

m 

DO  600  IV*1,NV 

ROLL 

122 

38 


FROG****  FINCON 


7*/75 


175 


ISO 


1*5 


150 


155 


200 


205 


210 


215 


220 


225 


1  FtNCON  7 +/7k  OPT«0  ROUND* •/  TRACE  FTN  4. A* Ail 

11/25/75 

11.(5. 

• 

ROLL 

IM 

•BEGIN  LOOP  OVER  SHIP  HEADING  (NUI,  PREDOMINANT  NAVE  DIRECTION. 

roll 

125 

• 

ROLL 

125 

00  500  INU*ltNNU 

ROLL 

125 

IF  (dV.EQ.l!  .AND,  (IMH.EQ.1H 

ROLL 

122 

2  NURCINU)  «  PI/100.  •  FLOAT  IHUOCIHUn 

ROLL 

12* 

OAHPUIINUI  *  OUCHV.il 

.OIL 

125 

IF  (ITEPATE  .EQ.OI  GO  TO  90* 

ROLL 

101 

NTRV  «  0 

.OIL 

lit 

yp  »  o.o 

.OLL 

102 

a  a  NTRt  •  nut  ♦  1 

.OLL 

101 

x  *  tp 

.OLL 

105 

M5FUIIMUI  •  OUCIIV.il  *  1.51*OUCUV.2>»X**(.772  ♦  l.M'DUGI M.S  ROLL 

105 

2)  •*  ♦  4.O0OUCIIV.4J*X»*2.0  ♦  9.40»DUCdVt5)PX”3.l  ♦ 

.OLL 

105 

2  24 • 0*DUC(IV » AI#X#*4.I 

.OLL 

102 

90  SIGLC(IHU)  *  SIGSOSCCIftJI  «  0. 

.OLL 

101 

SIGVLCdNJI  *  SGVSOSCdHUI  «  0. 

.OLL 

105 

• 

.OLL 

151 

•BEGIN  LOOP  OVER  SPREADING  ANGLE  <NUt. 

.OLL 

151 

• 

.OLL 

152 

•(FOR  PURELY  LONG  CRES TEO-CASE.  NO.  OF  NU'S  *  1.  ANO  NU  *  1U.I 

.OLL 

153 

4 

.OLL 

155 

IF  (ISC  •£ J.  01  NNU  *  1 

.OLL 

155 

00  400  INU*1.NNU 

.OLL 

155 

NUR(INU)  *HUR IIHU)4P(INUI  *COH 2 

.OLL 

152 

SINNU(INU)  •  SlNINURdNUII 

.OLL 

195 

50  COSHU(INU»  >  COSINUR (INUI  t 

.OLL 

155 

• 

.OLL 

Ml 

•BEGIN  LOOP  OVER  4AVE  FREQUENCY. 

.OLL 

2(1 

• 

.OLL 

M2 

00  200  IN*1.NW 

.OLL 

213 

KL  *  L  •  M( INI •  *2  •  COSNUdNUI  /  C2.  •  GRAVITY! 

ROLL 

2(5 

SINKL  *  SIR ( KL 1 

ROLL 

2(5 

COSKI  *  COS<KU 

ROLL 

2(5 

NE( IN* INUI * A9S ( M( INI  *  (1.  -  VC  1*0  •  VFSIIVI  /  GRAVITY  • 

ROLL 

2(2 

2  COSNUdNUMI 

ROLL 

((( 

TURF  •  HE (CM  « I NUI /MPHI 

ROLL 

((( 

BA  *  ? * • DAiPU (INUI  •  TUNF 

ROLL 

21( 

A  *  1.  -  TJNFMUNF 

ROLL 

211 

CA  *  SORT (A* A  ♦  BA*BAI 

ROLL 

(It 

• 

ROLL 

(13 

•TEST  FOR  MATERPLANE  SHAPE. 

ROLL 

(15 

* 

ROLL 

lit 

IF  (SHAPE  .EQ.  ELLIPI  GO  TO  110 

ROLL 

(15 

IF  (SHAPE  .EO.  RECTANG)  GO  TO  120 

ROLL 

(12 

* 

ROLL 

(11 

•MATERPLAVE  IS  PARABOLIC. 

ROLL 

(15 

• 

ROLL 

ttl 

LPRINE  *  L  *  .5 

ROLL 

221 

IF  ( IBS (KLt  .GT.  EPSI  GO  TO  105 

ROLL 

((( 

F  *  1. 

ROLL 

223 

G  *  0. 

ROLL 

2(5 

0*1. 

ROLL 

ttt 

GO  TO  127 

ROLL 

(15 

105  F  *  3.  •  (SINKL  -  KL  •  COSKLI  /  KL**3 

ROLL 

(12 

G  *  1575.  /  KL**7  •  (d.  -  2.  •  KL**2  /  5.!  •  SINKL 

ROLL 

ttl 

2  -  KL  •  (1.  -  KL"2  /  15.1  •  COSKLI  -  F 

ROLL 

(tt 
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PROGRAM  FIMCON 


7 k/7k  OPT»0  ROUND**/  TRACE 


U/H/M  ll.7t.lt 


710 


71* 


7*0 


7** 


7*0 


75* 


700 


76* 


770 


77* 


700 


70* 


0«F*MN/GN«G 

ROIL 

71* 

GO  TO  127 

ROLL 

711 

# 

ROLL 

717 

•MATERPLANf  IS  EL.IPTICAL.  {DETERMINE  THE  ZEROTH  AND  FIRST  MESSEL 

ROLL 

711 

•FUNCTIONS  OF  KL.I 

ROLL 

716 

• 

ROLL 

71* 

110  t PRIME  •  SORT  1 7. 1/* •  •  L 

ROLL 

716 

IF  IIBSIKLI  . GT  •  EPS!  GO  TO  115 

ROLL 

297 

F  «  1. 

ROLL 

238 

G  »  0. 

ROLL 

239 

0*1. 

ROLL 

290 

GO  TO  177 

ROLL 

291 

It*  CALL  NESSJ  IXL.0.,1, BESSEL! 

ROLL 

292 

F  «  7.  /  XL  •  BESSEL  171 

ROLL 

293 

G  *  -  ».  /  KL**7  •  BESSEL  111  *  10.  /  KL**7  '  1.1  •  F 

ROLL 

299 

0»F40H/GH»G 

ROLL 

76* 

GO  TO  177 

ROLL 

296 

9 

ROLL 

767 

ROLL 

298 

• 

ROLL 

769 

120  1* PRIME  «  t 

ROLL 

75* 

IF  IA5SIKLI  .GT.  EPS!  GO  TO  17* 

ROLL 

7*1 

0*1. 

ROLL 

7*7 

GO  TO  127 

ROLL 

293 

125  0  *  SINKL  /  KL 

ROLL 

299 

177  CONTINUE 

ROLL 

7*5 

KLPRINE  *  .5  •  N(IMI**7  /  GRAVITY  •  L PRIME  •  tOSNUlINUI 

ROLL 

7*6 

C  *  SINIKLPRINEI  /  XLPRIME 

ROLL 

7*7 

H  «  0  -  Q*C*TUNF*TUNF 

ROLL 

7*6 

9 

ROLL 

75* 

•COMPUTE  ROLL  TRANSFER  FUNCTION. 

ROLL 

76* 

♦ 

ROLL 

761 

TRIM!  *  E<TIIMI/CA*SIHNU(INU1*S0RTIH*H*C*C»BA*BAI 

ROLL 

767 

* 

ROLL 

761 

•COMPUTE  ROLL  SPECTRUM. 

ROLL 

766 

♦ 

ROLL 

76* 

SURIIM.INUI*  SIIMI»TRIIMI»TRIINI 

ROLL 

76* 

SURVIIH.INUI  •  SURIIN,INUI*HE(IN>INUI>NE(IN.INUI 

ROLL 

767 

9 

ROLL 

76* 

VEND  OF  LOOP  OVER  WAVE  FREQUENCY, 

ROLL 

76* 

• 

ROLL 

77* 

200  CONTINUE 

ROLL 

771 

• 

ROLL 

777 

•DETERMINE  R*S  ROLL  VALUE. 

ROLL 

771 

• 

ROLL 

776 

CALL  ALGRNS  INN.N.  SUR  (1.  INUI  .  MOUR 1 1 NUI 1 

ROLL 

77* 

CALL  ALGRNG I NN.M . SURV (1 . INUI  .VMOURIINUII 

ROLL 

776 

♦ 

ROLL 

777 

•  IF  REOUESTEO  compute  lomg-cresteo  VALUES 

ROLL 

77* 

• 

ROLL 

77* 

IF  use  .ME.  01  GO  TO  310 

ROLL 

7*1 

SIGLCIIMUI  *  SQRTIMQURtlNUII 

ROLL 

7*1 

SIGVLCIIMJ)  *  SORTI VMOUR 1 IMJ 1 1 

ROLL 

7*7 

9 

ROLL 

7*1 

•  IF  REQUESTED,  ITERATE  OVER  ROLL  DAMPING  FOR  LONG-CREST  CASE. 

ROLL 

7*6 

• 

ROLL 

7** 

IF  1  ITERATE  .EQ.  01  GO  TO  600 

ROLL 

7*6 

40 


PROGRIN  FI  (ICON 


74/7*  OPT»l  ROUND**/  TRICE 


FT#  4.6*461 


11/M/M  11.21.11 


KHJ  »  INU 

ROLL 

207 

KV  »  IV 

ROLL 

211 

KNU  «  INU 

ROLL 

211 

S6V0MLC  «  SIGVLCI IHU)/*PMI 

ROLL 

211 

2H 

PHIN  *  04N»U<INUI 

ROLL 

211 

CALL  ITREQ(SG VOULC9  *  RETURNS!  OOoMliOOO) 

ROLL 

212 

• 

ROLL 

213 

•KG  IN 

SUNNING  Of  SHORTCRCSTEO  RESPONSE  OATA. 

ROLL 

214 

• 

ROLL 

211 

211 

390 

SIG SQSC ( INJI  *  SIG$0$C(INU>  ♦  COS1 (INU)*H|U* (INU 1 

ROLL 

216 

SGVSCSC (INU)  »  SGVSQSC (IHU)  ♦  COS 1 (INU) 9VN|UR (INU) 

ROLL 

217 

• 

ROLL 

211 

UNO  Of 

■  LOOP  OVER  NU. 

ROLL 

211 

• 

ROLL 

311 

111 

Ml 

CONTINUE 

ROLL 

311 

• 

ROLL 

302 

IF  (ISC  .EQ.  II  INU  *  1 

ROLL 

313 

If  (ISC  *E3.  0)  GO  TO  210 

ROLL 

314 

SIGSC(INU)  a  SORT  (CONlMIGSQSC  ( IHU)  1 

ROLL 

311 

111 

St GVSC(IHJ)  »  SORT (CONI* SGVSQSC (INU) 1 

ROLL 

316 

SGVONPH  «  SIGVSC (INUI/MPHI 

ROLL 

317 

RNU « I HU 

ROLL 

311 

KNU  « INU 

ROLL 

311 

KVa  IV 

ROLL 

311 

111 

If  (ITERATE  .EQ.  0)  GO  TO  250 

ROLL 

311 

PNlNaOAMPU(INU) 

ROLL 

312 

CALL  ITRE)(SGVONPHI'  RETURNS (00*2509500) 

ROLL 

313 

* 

ROLL 

314 

•IF  REQUES ICO .  PRINT  LONGCRESTEO  SPECTRUM  IND  ITS  COHPONCNTS. 

ROLL 

311 

111 

• 

ROLL 

316 

211 

IF  (IPRINTIll.NE.SPCCTRIt  GO  TO  211 

ROLL 

317 

• 

ROLL 

311 

IF  ( IUN1TS  .EQ.  II  GO  TO  211 

ROLL 

311 

• 

ROLL 

321 

321 

NRITE  (6,30111  T1TLE2,SHHETIIHHI,TIIImHI.«RII«I 

ROLL 

321 

GO  TO  232 

ROLL 

322 

• 

ROLL 

323 

231 

NRITE  (6,20111  TITLE2.SNN  IIHH 1 .  THINK  1  *  VKIZVI 

ROLL 

324 

232 

CONTINUE 

ROLL 

321 

321 

NRITE  (6,  21231  NUDIINUI 

ROLL 

326 

NRITE  (6,21201 

ROLL 

327 

DO  310  IN*1,NH 

ROLL 

321 

PER  •  2.»PI/N(INI 

ROLL 

321 

LAN  »  PER*»ER*GR6WY/(Z.»PII 

ROLL 

331 

331 

LINONL  >  LAN/L 

ROLL 

331 

NN  •  360. *4 (INI *N(INI/(2.*PI*6RIVITTI 

ROLL 

332 

rao*tr(ini»tr(ini 

ROLL 

333 

SO  •  S ( INI/ ( NN*NNI 

ROLL 

334 

RIOO  »  NN*NN*RIO 

ROLL 

331 

331 

SIRO  •  SO  •  RIOO 

ROLL 

336 

NRITE  (6,21211  N(INI  .NEdN.INUI  .LIN.LINONL.NN.TRdNI  .RIO.  St  INI* 

ROLL 

337 

2 

SURTINI .  RIOO.SO.SURO 

ROLL 

331 

350 

CONTINUE 

ROLL 

331 

NRITE  (6,21221  SIGLCdNUI 

ROLL 

341 

341 

• 

ROLL 

341 

•IF  REQUES TEO,  COMPUTE  STISILIZED  ROLL  IND  FIN  NOTIOAS  FOR  THIS 

ROLL 

142 

•predominant  heroins. 

ROLL 

343 

41 


PBOGBAH  FINCON 


74/71.  OPT«0  ROUNO* •  /  TRACE 


FTN  4.6*660 


11/20/79 


345 


350 


355 


360 


365 


370 


375 


300 


385 


390 


393 


250  IF  INSTAB  .LT.  II  GO  TO  500 
DO  280  IS*1.NSTA» 

CALI.  FINSTAB  (151,  RETURNS  (280,9001 
280  CONTINUE 

• 

TEND  OF  LOOP  OVER  PREOOHINANT  HEADING, 

500  CONTINUE 

• 

'ONLY  LONGCRESTEO  ROLL  VALUES  ARE  OUTPUT, 
v 

IF  (IUNITS  ,EO.  01  GO  TO  502 

WRITE  (6,39101  TITLE2 .SWHET (I WHI • TO (IWHI «  VK  ( IVI 

GO  TO  503 

• 

502  WRIT£(6, 20101  TI TLE2, SMH ( IWHI  , TO (I WHI »VK (IVI 

503  CONTINUE 

WRITE  111  TITLE 2 , S WHI IWHI , TO (IWHI ,VK (IWHI 

IF  (ISC  ,N£.  01  GO  TO  550 

IF  (ITERATE  • EQ«  01  WRITE (6,2 0301 

HRITEI6, 20111 

(10  505  I  HU  •  1,NHU 

505  WRITEI6. 20141  HUO I IHUI, OAHPUI  IHUI , S IGLC (IHUI 
WRIT'lll  (51 GLC  1 1  HU  I  ,  INU  •  1,  NMUI 
IF  INSTAB  ,EQ •  01  GO  TO  600 
00  540  I  •  I.NSTAB 
WRITEI6. 20321  I 
00  525  I  HU  «  l.NHU 

525  HRITEI6, 20141  HUO I IHUI ,0AHPSI I • I HU I ,SS IGLCI I , IHUI ,1 TECT ( I, IHU I , 
2  BMTlCII. IHUI  .BVELLCII, IHUI 

WRITE (11  (SSI GLC ( I, IHUI •  IHU  «  l.NHUI 
540  CONTINUE 

IF  (NSAT  .£9,  01  GO  TO  600 
WRITEI6, 22721 
GO  TO  600 

m 

•ONLY  SH0RTCRESTE3  ROLL  VALUES  ARE  OUTPUT 
• 

550  IFI  ITERATE  ,EO.  01  WRITEI6, 20301 
WRITE  1 6, 20111 
00  555  IHU  ■  l.NHU 

555  WRITEI6, 20141  HUD  I IHUI .OAHPUI I NUI .SIGSCIIHUI 
WBITEdl  (SIGSCIIHUI,  IHU  •  l.NHUI 
IF  INSTAB  ,EQ.  01  GO  TO  600 
00  590  I  >  I.NSTAB 
WRITEI6, 20261  I 
00  575  IHU  ■  l.NHU 

575  HRITEI6, 20141  HUO (IHUI .OAHPSI I, IHUI .SSIGSCI I. IHUI . ITEST (I, IHUI . 
2  BHOT SC  1 1  * INU I , BVELSC 1 1, IHUI 

WRITE  111  ISSIGSCII, IHUI, IHU  >  l.NHUI 
590  CONTINUE 

IF  (NSAT  .19.  01  GO  TO  608 
WRITEI6, 22221 
F 

SEND  OF  LOOP  OVER  SPEEO. 


ROIL 

346 

ROLL 

345 

ROLL 

366 

ROLL 

347 

ROLL 

340 

ROLL 

349 

ROLL 

350 

ROLL 

351 

ROLL 

352 

ROLL 

353 

ROLL 

356 

ROLL 

355 

ROLL 

356 

ROLL 

357 

ROLL 

350 

ROLL 

359 

ROLL 

360 

ROLL 

361 

ROLL 

362 

ROLL 

363 

ROLL 

366 

ROLL 

365 

ROLL 

366 

ROLL 

367 

ROLL 

360 

ROLL 

369 

ROLL 

370 

ROLL 

371 

ROLL 

372 

ROLL 

373 

ROLL 

376 

ROLL 

375 

ROLL 

376 

ROLL 

377 

ROLL 

370 

ROLL 

379 

ROLL 

300 

ROLL 

301 

ROLL 

302 

ROLL 

303 

ROLL 

306 

ROLL 

305 

ROLL 

306 

ROLL 

307 

ROLL 

300 

ROLL 

309 

ROLL 

390 

ROLL 

391 

ROLL 

392 

ROLL 

393 

ROLL 

396 

ROLL 

395 

ROLL 

396 

ROLL 

397 

ROLL 

390 

ROLL 

399 

ROLL 

400 

42 


PROGRAM  FIN  CON 


76/7*  DPT»0  SOUND-  •/  TRACE 


FTN  *.6**61 


11/20/79  a.tf.ii 


Ml 

% 

ROLL 

*01 

Ml 

CONTINUE 

ROLL 

*02 

• 

ROLL 

*03 

9ENO  OF  LOOF  OVER  SE*  CONDITION. 

ROLL 

*0* 

• 

ROLL 

*09 

All 

700 

CONTINUE 

ROLL 

*06 

900 

HRITE  (6.20191 

ROLL 

*07 

• 

ROLL 

*00 

STOP 

ROLL 

*09 

1000 

FORMAT  (0A10I 

ROLL 

*10 

*11 

1001 

FORMAT  (16191 

ROLL 

*11 

1302 

FORMAT  IAF10.6I 

ROLL 

*12 

lwOt 

FORMAT  (I9,2F10.9,A10» 

ROLL 

*13 

2000 

FORMAT  1  INI.  27 ( /I **9X.*2M*  •  •  ROLL  NOTION  PREDICTION  PROGRAM  •  • 

ROLL 

*1* 

2 

•Z//.96X, 3*101 

ROLL 

*19 

*1* 

2001 

FORMAT  (1H1.GA16///! 

ROLL 

*16 

2002 

FORMAT  (IX,* SIGNIFICANT  NAVE  MCICMT (SI  (FEET!  »*,9X,  10FT.2I 

ROLL 

*17 

2003 

FORMAT  (IX, •  MODAL  NAVE  PERIOO(S>  (SECONDS!  -*.*X,1*F7.2//1 

ROLL 

*10 

2000 

FORMAT(/1X,*OISPLACEMENT  (L.  TONS!  **16X,F7.0/1X  *LEMGTM  SETREEN 

ROLL 

*19 

2PP  (FEETI  ••1*3. FT . 1/ IX. -OR AFT  (FEET!  ■•  26X.F7.2/1X, 

ROLL 

*20 

»n 

2- TRANSVERSE  NET ACENTRIC  HEIGHT  (FEETI  ••2X.FT. 1/1X*NETACENTER  A NOV 

ROLL 

*21 

2E 

:  BUOTANCV  CENTER  (FT!  «• . 1X.F7 . 2/lX,*ROLL  PERIOD  (SECONOSI  «*,1TX 

ROLL 

*22 

2. FT.  2/it.  *9  **,37X,F7.3I 

ROLL 

*23 

2009 

FORMAT  (1X,*NATEAPLANC  SMAPE  •*,2*X,A10I 

ROLL 

*2* 

2006 

FORMAT  (IX, ^SPREADING  ANGLE  ■•23X.F7.0//! 

ROLL 

*29 

*21 

2007 

FORMA  T ( IX.  110MDAMPING  INPUT  IN  THE  FORM  N  »  Cl  ♦  1.61*C2*T**0.7 

ROLL 

*26 

272  ♦  l.«0*C3*V  ♦  *.00*C**V**2  ♦  9.*t*CS*T**3  ♦  2*. 30*C6*T***// 

ROLL 

*27 

2" 

1  SPcEO  (KNOTS »*.*X.*C1*.TX.*C2*,TX,*C3*.TX,*C**.TX.*CS*,TX.*C**I 

ROLL 

*20 

2009 

FORMAT  (1X,*SMIP  SPEEDISI  (KNOTS!  •*.  19X.9F9 . 1/1 

ROLL 

*29 

2009 

FORMAT  (*X,  F9.1,7X,6(2X«F7.*I  1 

ROLL 

*00 

*3« 

2010 

FORMAT  (1M1.SA10//1X.MS1GH1FICANT  HAVE  HEIGHT  •»F7.2.*  FEET* 

ROLL 

*31 

2 

/  IX  *NOOA 

ROLL 

*32 

2L  NAVE  PERtJD  ••FT.**  SECONDS- /IX*  SMI*  SPEED  ••  FS.l,*  KNOTS*  1 

ROLL 

*33 

2011 

FORMAT  (//1X,*UNSTA0ILI2ED  RMS  ROLL  (DEGREESI*/1X.*MEA0ING*7X*LC* 

ROLL 

*3* 

2 

OX, *SC*/I 

ROLL 

*39 

*39 

2012 

FORMAT  (I6.1X,2F1*.2.2X,*1,19X,2F1*.2.12X,2F10.2! 

ROLL 

*36 

2013 

FORMAT  (//IX » *UNSTABILI7EO  RMS  ROLL  (OECREESI*/  IX, -HEADING*  OX. 

ROLL 

*37 

2 

•N*8X*SC*/I 

ROLL 

*30 

201* 

FORMAT  (I6.1X,F10.*,F10.2,2X,*1,29X.F10.2,22X,F10.2I 

ROLL 

*39 

2019 

FORMAT  I1M1,27(/I  ,97X,17H*  **ENO***! 

ROLL 

**■ 

**0 

20  1  E 

FORMAT  (IX,  *ROLL  DAMPING  COEFFICIENT(S!  •*,  13X ,SFS. 3///I 

ROLL 

**1 

2017 

FORMAT  ( IX, -  ITERATION  OVER  ROLL  DAMPING  HILL  OE  OONE.VI 

ROLL 

**2 

2010 

FORMAT  (//1X,*L0NGCRESTED  SPECTRA  AND  COMPONENTS  MILL  IE  PRINTEO. 

ROLL 

**3 

2M 

ROLL 

*•* 

2019 

FORMAT  ( //IX ,* INTERMEDIATE  STEPS  IN  ROLL  DAMPING  ITERATION  HILL  • 

ROLL 

**! 

**5 

2t 

:  PRINTEO. *1 

ROLL 

**6 

20  20 

FORMAT  (//9X,*N*,0X,*NE*,7X,*LAM*,9X,  *LAN/L*,9X,*K*,0X, 

ROLL 

**7 

2 

•TR*,7X,*RA0*,*X,»M  SL  S*,7X,*SUR*,  3X,*RAO  OIM*,*X,*M  NT  S*,3X, 

ROLL 

*** 

2 

•SUR  DIM*/  3X,»RA0/SEC*3X  *RAO/SEC*OX.*FT* . 17X,*OEG*,3X,*OCC/OEG 

ROLL 

**9 

2 

DES/DEG  so  DEGSQ  SEC*1X.*DEGS0  SEC*.1X,*0EG/FT  SQ*1X,*FTS0  $EC*1K 

ROLL 

*90 

*90 

2 

•OEGSQ  SEC*! 

ROLL 

*91 

2021 

FORMAT  (12F1(,]| 

ROLL 

*92 

2022 

FORMAT  ( //IX , -RMS  ROLL  »*F7.2*  DEGREES*! 

ROLL 

*93 

2023 

FORMAT  ( //IX , -SHIP  HEADING  **I3*  DEGREES*) 

ROLL 

*9* 

2026 

FORMAT  (//1X,*R0LL  STABILIZATION  MILL  9E  CALCULATED  FOR  *13 

ROLL 

*99 

*99 

2 

*  CASES*//  IX  *FIN  AND  CONTROL  STSTEN  PARAMETERS  ARE  AS  FOLLONSI* 

ROLL 

*90 

2 

//•  CASE*  2X  *M*  6X  •*»  5X  *R  IDCL/OBIFS*  *X  *M0*  *X  *M1*  *X 

ROLL 

*97 

43 


PROG* AH  FINCOM 


74/7*  OPT'O  90UN0"/  TRACE 


FTN  4.6*660 


tl/n/M  ll.25.JI 


460 


465 


470 


675 


400 


605 


490 


7  'M2'  40  »HT»  60  'M4'  4*  'GK»  4*  »GV*  4*  *K1*  60  4*  *KS*  40 

7  'At •  4*  »A2'  40  'A3'  40  '81*  40  40  *63*  /10K  *PT  SO*  40  *FT* 

7  40  'PER  OE5'/) 

7025  FORMAT  (15.I3.F7.2.F6.2.F12.3.16F6.3) 

2076  FORMAT  < //IX , 'C ASE'I 3 .» (  STA9ILIZC0  RMS  ROLL  (OEGRCCSI* 

2  100.*FIN  MOTION  (DEGREES)*  tOX  *FIN  VELOCITY  I DEGREE S/SECOND) •/ 

2  10* •HEADING*  AO  *N*  00  »SC*  360*SC*  300  *SC*/» 

2027  FORMAT  I//1 X . *CAS£*I3,*I  STABILIZED  RMS  ROLL  (DEGREES)* 

2  100  ,'FI6  MOTION  (DEGREES)*  ISO  *FIN  VELOCITY  (DEGREE S/SECOMO)  •/ 

2  10*  *Mf  AOINS*  70  *LC*  AO  'SC'  270  *LC*  50  *SC*  190  *LC*  00  'SC*/) 

7020  FORMAT  (//5X**IV*50, *BST0P*5  0*  *BVELMAO*/* 

2  (I7,7X,F10.5.20,F10.5,/) ) 

2029  FORMAT <//10**MO  SPECTRA  OUTPUT  AILONED  IN  SMORTCRESTEO  CASE.  REV 
21 SE  YOUR  IN»UT  DATA.*) 

2030  FORMAT(//10,«ITERATION  OVER  ROLL-OANPIMG  NOT  PERFORMED.*) 

2031  FORMAT  (//10,*UMST  ABILIZEO  RMS  (XL  (DEGREES) */10 ,*MEA01NG*, 00  . 

2  *N*,00,*LC*/( 

2032  FORMAT ( //10 ** CASE *13**)  STABILIZED  RMS  ROLL  (DEGREES)* 

2  100 , 'FI N  MOTION  (OEGREESI  *»100,*FIN  VELOCITY  (DEGREES/SECOND) */ 

2  io,*meaoin;*  ax  *n*  ox  *lc*  J60  *LC*  SOX  *LC*/) 

2222  FORMAT  (//10.36H*  SATURATION  FACTOR  IS  INSIGNIFICANT) 

3002  FORMAT  (10, 'SIGNIFICANT  NAVE  HEIGHT (S)  (METERS)  *',SO,10F7.2) 

3004  FORMAT(/1k,'OISPLACEMENT  (N.  TONS)  "230, FT.  0 /IX  ‘LENGTH  BETWEEN 
2PP  (METERS)  **140, F7. 1/10. *ORAFT  (METERS)  »•  260, F7. 2/10, 

2* TRANSVERSE  METACENTRIC  HEIGHT  (METERS)  "TX.FT.l/lO'HFTACENTE*  OB 
20VE  BUOYANCY  CENTER  (M)  *•  ,40. F7. 2/10, 'ROLL  PERIOD  (SECONOS) 

2190, FT. 2/10*  *0  "*39X,F7.3> 

3010  FORMAT  (1H1.0A10//1X, 'SIGNIFICANT  NAVE  HEIGHT  "FT. 2,'  METERS' 

2/10  'HOOA 

2L  NAVE  PERI30  "FT. 2'  SECONDS' /1X*SHIP  SPEED  »•  F5.1,*  KNOTS') 

3024  FORMAT  (//IX. 'ROLL  STABILIZATION  MILL  9E  CALCULATED  FOR  *13 

2  •  CASES'//  10  'FIN  AND  CONTROL  SYSTEM  PARAMETERS  ARE  AS  FOLLONSI' 
2  //•  CASE'  20  'M'  40  'A*  70  »R  (DCL/OBIFS*  40  'HO'  40  'Ml*  *0 
2  'M2'  40  'M3'  40  'M4'  60  'GK'  4X  'GV'  40  'Kl«  40  •«'  40  'KS'  40 
2  'Al'  40  'A2'  4X  'A3'  4X  'BI'  4X  '92'  40  'B3»  /110  'M  SO'  5X  'H  * 

2  3X  'PER  OEJ'/) 

END 


ROLL 

45» 

ROLL 

459 

ROLL 

460 

ROLL 

461 

ROLL 

462 

ROLL 

463 

ROLL 

466 

ROLL 

665 

ROLL 

666 

ROLL 

66T 

ROLL 

666 

ROLL 

669 

ROLL 

4TB 

ROLL 

4T1 

ROLL 

*72 

ROLL 

473 

ROLL 

6T* 

ROLL 

4TS 

ROLL 

476 

ROLL 

677 

ROLL 

47B 

ROLL 

479 

ROLL 

*00 

ROLL 

461 

ROLL 

662 

ROLL 

663 

ROLL 

666 

ROLL 

665 

ROLL 

666 

ROLL 

667 

ROLL 

666 

ROLL 

469 

ROLL 

690 

ROLL 

491 

ROLL 

492 

ROLL 

493 

ROLL 

496 

44 


SUBROUTINE  ITREQ 


FTN  t.i»UI 


U/2A/7S  U.M.K 


mmmmm 


IS 


IS 


20 


2S 


31 


74/71.  OPT*0  NOUNO* •/  TRACE 


SUBROUTINE  I TREO  1SI6I.RE TURNS l»,S.C» 

CONNON/mTN/STAT,PRCN,NTRV,  I NU, INU.IN.CR.CRNl.TRNl.TR, VPP1, 
PHIN,  IPRINT  I  SI,  ITERATE, NPHI 
OATS  ITERATN/7HI TERATN/ 

IF  INTRV.3T .111  NRITEIA.440SI 
IF  (NTRT  .ST.  HI  RETURN  C 

.eq.iteratni  .and.  iinu  .eq.  ii  .ano.  inert. eo. hi 


•  EO.  ITERATNI  NRITEIA.SISSIIV.INU.NTRV.PNIN.VPPI, 


(CP  -  GPNlt I 


IF  (  (IPRINTI2I 
2  HRITEIA, 11011 
EP  >  STAT*SIG 
IF  ( I PRINT  121 
2  VP,VPN1,G».GPN1 

IF  ISP  .LT.  .All  RETURN  S 
IF  INTRV  .EO.  11  SO  TO  20 
IF  IASSI1  -  GP/TPI  .LE.  PRCNI  RETURN  3 
TPP1  «  (G>N1*TP  -  GP*YPNll/t ITP  -  TPN11 
t»ni  *  vp 
GPHl  •  CP 
VP  »  VPP1 

IF  I IPRINT 1 2  I  .EO.  ITERATNI  NRITEtt.SlIOlIV.INU.NTRV.PNIN.VPPl, 

2  VP  .VPN1.6P.6PN1 
RETURN  A 
20  GPN1  «  GP 

VPN1  «  0.1 
VP  •  GP/2.0 

IF  I IPRINT 121  .EO.  ITERATNI  NRITE IS. 01001 IV.INU.NTRV.PNIN. VPPl, 

2  VP,VPN1,G»,GPN1 
RETURN  A 

4S98  FORNAX  I  IN  ------  PROGRAN  STOPPED  BECAUSE  ROLL  FAILED  TO  CONS 

2ERGE  NITMIN  10  ITERATIONS  -  -  -  -  -  M 

5100  FORNATI1X, SIS, 6F 10.31 

5101  F0RNATI//,SX,»I»«,3X,*INU*.2X,»NTRV»,%*.PPII1NP,AX,»VP»IP,7»,»VP» 
2,7x,*VP-l«.7X,*GP*,7X,»GP-l*/l 

ENO 


TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 

TREO 


oo  u  oo  u u  og  gg 


SUBROUTINE  BRUSSP 


7*/7*  orr«»  rounob*/  trace 


FTN  T.IMtl 


11/II/M  ll.tf.lt 


1 


If 


SUBROUTINE  fRMSSP  TN.SISNM. TNOOAL, N.St  BBSS 

■isitstmtuiiamiinuimaimmui^niunaHnuinimiunHii  BNSS 
COMPUTES  BRET SCHNEIDER  NAVE  SLOPE  SPECTRUM  BBSS 

REEL  K  BBSS 

OIHENSION  4<N I , SINt  BNSS 

OAT  A  A.B.PI.G  /tlT.KNtlfW.tfHil.UlMirilMm/  BBSS 

«  TNOOAl  >  2.7*  •  SIGHM**.S  t  BBSS 

TNOOAL*  •  TNOOAL***  BNSS 

. BBSS 

FOR  PIEPSON»NOSKONITZ  SPECTRA  -  TNOOAL***  ■  ft. IBS*  *  SICBH**2  BBSS 
CONI  *  A  •  SI6NM**2  /  TNOOAL*  BNSS 

CON  2  *  8  /  TNOOAL*  BNSS 

00  It  IbI.N  BNSS 

N*  b  Mill***  BBSS 

N5  «  NUT  •  N*  BNSS 

ARS  •  COM2  /  N*  BBSS 

if  <arc  .or.  sio.i  sat  b  «.  bbss 

IF  (ARC  .ST.  5 to. I  SO  TO  It  BNSS 

■*•■*■******•■•■*•■»■*****■■■■■  ■•■■•••••MBBB*«BBM««BSB*«BBBB»BBMBB««ia  BBSS 


BRET SCHNEIDER  BABE  HEIGHT  SPECTRUM 
Slit  •  CONt/NS  •  EXPI-ARGI 

lasaaaaaa aaaBBaaBBBaa BBaaaaaaaa  asaaaaaaaaaaBBas 

HAVE  NUMBER  IN  DEGREES 
K  >  Itt.*H(II*Hat  /  T2.*PI*G  I 


ORE T SCHNEIDER  HAVE  SLOPE  SPECTRUM 

siit  b  k*k  •  sat 

It  CONTINUE 
RETURN 
ENO 


BBSS 

BBSS 

laaaaaaaaBaaaaaaaaaaaaa  BBSS 
BBSS 
BBSS 

rBBBBBBBBBBBBBBBBBBBUB  BBSS 

BBSS 

BNSS 

BBSS 

BBSS 

BBSS 
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subroutine  ncrnc 


71. /7k  OPT»0  POUND**/  TROCE 


FTN  k.(«kU 


11/M/M  11.25.30 


t 

c 

c-« 

- VERSION  3  -  COC  0700  -  0L6RNG  - 

JONUORT.  1970 - 

OLGR 

OLGR 

2 

S 

c 

OLGR 

5 

* 

SUBROUTINE  OLGRNC  (N.N.S.OREOI 

OLGR 

6 

♦ 

ALGA 

f 

•this  subroutine  computes  the  oreo  unoer  the  curve 

FOR  0  PORTICULOR 

OLGR 

• 

•SPECTRUM*  IN  ODO  NUMBER  Of  POINTS  (FREQUENCIES) 

SHOULD  BE  USED.— — — 

OLGR 

5 

• 

OLGR 

II 

10 

INTEGER  ERROR 

OLGR 

11 

DIMENSION  4(H)* SIN) 

OLGR 

12 

DATA  I T AG/D/ 

OLGR 

13 

ooro  EPS/O* 0  000000501/ 

OLGR 

it 

c 

OLGR 

15 

19 

ERROR  •  10 

OLGR 

to 

IF  (ERROR  *EQ •  101  ITRG  >  0 

OLGR 

ir 

f  RR3R*0 

OLGR 

to 

mo*vio-oO) 

ALCR 

15 

ARE  AO  »  0.5*  S  1 1  •  •  OMl  1  -MQI 

OLGR 

21 

70 

HN«N-2 

OLGR 

21 

AREA«0. 

OLGR 

22 

T£<fP  *  0. 

OLGR 

is 

ARE  A 2  *  AREAS  *  0. 

OLGR 

2N 

N ONEGA  «  100tN»2> 

ALtft 

25 

25 

DO  70  H.1.H5.2 

OLGR 

20 

AsN(«A2»-W('U 

OLGR 

27 

8«N<4a2»-NHA1) 

OLGR 

20 

OLGR 

25 

PURE*  *  0*0/0. *(SIHI*l3.*C*OI/IO*Ct*SI5MI*0/(B*CI* 

OLGR 

31 

30 

2  S(5*2l*(2.**-3.*Ct/(0*tl< 

OLGR 

31 

TENP  »  PAREA 

OLGR 

32 

IF  (PORE*  *LT.  0.1  IENP  •  0. 

OLGR 

33 

AREA  *  AREA  A  TENP 

OLGR 

30 

IF  (PORE*  *GE.  0.1  GO  TO  20 

OLGR 

35 

35 

I FI* PORE*. GT. 0 *10* ORE At  ERROR* 1 

OLGR 

30 

20 

CONTINUE 

OLGR 

37 

PORE*  *  0. 

OLGR 

30 

ir  (ORE*  ..E.  0.000010)  GO  TO  100 

OLGR 

35 

c 

OLGR 

NO 

Itlllv"  '  W  Srtli  IKAl  0 LU)U Kt*” 

c 

OLGR 

hi 

TEMP  «  0. 

Al»6t 

NS 

SNA t  *  S(l> 

OLGR 

ITEST  «  1 

OLGR 

N5 

*5 

DO  30  I»2.N 

OLGR 

hi 

30  IF  (SHI  .ST.  SHOXI  SMX  «  Sill 

OLGR 

ST 

X10PRCN  *  O.IO'SNOX 

OLGR 

NO 

50  ITEST  »  ITEST  *1 

OLGR 

N5 

IF  (ITEST  .EO.  2)  J*1 

OLGR 

50 

IF  1  ITEST  .EO.  31  J*N 

OLGR 

51 

IF  ((SM0X*SIJ)l  .LE.  EPS  1  ERROR  *  ITEST 

OLGR 

52 

IF( (ERROR  *EO.  21  .OR.  (ERROR  .E0.3II  ITOG 

■  1 

OLGR 

53 

IF! (ERROR  .EO.  2)  .ONO.  1 ITEST  .EO.  2)1  TEMP* ORE O.OREOO 

OLGR 

SN 

IFI (ERROR  .EO.  21  .ONO.  (ITEST  .EO.  211  OREO*TEHP 

OLGR 

55 

55 

IF  (S(JI  .JT.  X10PRCNI  ERROR  »  ITEST  ♦  2 

OLGR 

to 

IF  ((ERROR  .CT.  II  .ONO.  (ITOG  .EO.OII  ITOG 

■  1 

OLGR 

57 

IF  (IJ  .CO.  Nl  .ONO.  (ERROR. EO.  Nil  60  TO  117 

OLGR 

50 

47 


-- ri‘i<»  r*-'  • 


SUBROUTINE  ALGRNG  74/74  OPT«0  ROUND**/  TRACE 


FTN  4.6*461 


u/ii/ti  11. is. si 


IF  ((J  .EO.  NT  .ANO.  TERROR  .IT.  Sit  SO  TO  111 
IF  ((ERROR  .EQ.4I  .OR.  (ERROR  .CO.  SMSO  TO  II 
IF  (ITEST  «LT .  31  60  TO  SO 
60  TO  111 

-ORAM  A  STRAIGHT  LINE  THRU  FIRST  (LASTI  TMO  SPECTRAL  VALUES 
-TO  THE  ASSIISSA  ANO  AOO  ON  AREA  FOR  CLOSURE  AT  LON  (HI6HI 
-FREQUENCY  ENO .  — — - — - - - ............. 


60  IF  ((ERROR.EQ.SI  .ANO.  (S<JI  .SC.  SIJ-1MIG0  TO  7S 
IF  (IERR0R.Ea.4l  .ANO.  (Sdl.SC.  S(t  11160  TO  SO 
IF  (ERROR  .£0.  41  J*I 
IF  (ERROR  .EO.  St  J*N 

SLOPE  •  tS(J-lt  -  S(JII/(N(J-1I  -  MtJII 
IFI SLOPE  .LE.O.I  60  TO  TO 
SLOPE  •  AHIN1 (-SLOPE. -1.01 
70  IF  IJ  . EQ.  It  J  *  1 

INCH  •  -0.1  •  S(JI**I  /  SLOPE 
IF  (ERROR  .EQ.  41  AREAS  «  ANI N1 (ARCAI.ANEMI 
IF  (ERROR  .EQ.  SI  AREAS  *  ANEM 
7S  TENP  *  AREA  ♦  AREA!  •  AREAS 
IF  (J  .LT.  Ml  60  TO  SI 
117  TEMP  *  AREA  ♦  AREA!  ♦  AREAS 
AREA  «  TEN* 

111  RETURN 


SUBROUTINE  FIRSTAR 


FTN  k.6»ktO 


II /Zt/M  II 


f ,  vK%.  -  , 


7A/7A  O*T*0  ROUND**/  TRACE 


1 

SUBROUTINE  FINSTAB  IISI,  RETURNS  UAA.Mtl 

FST 

• 

FST 

•SUBROUTINE  TO  PRiOICT  STHUIIEO  (FINI  ROLL  NOTION  USING  CONOLLT  AND 

FST 

•CO*. 

FIT 

5 

• 

FST 

REAL  NN.K<U,««Z,KK3,KR,(I,K1,KZ,K3 

FST 

CONNOR/ ITRTN/STAT  «PRCN,NTRV« INU.KNU, I T, GF.GPNl , TPN1 ,TP, T..I, 

FST 

2 

PMIR.IPBIRTH)  , ITERATE, NPHI 

FST 

CONNOR  /STAB/  NSTAB.HtlQt  ,  AREA  CIO  1  «Rt 101  ,OCLO tFS till  «Hlf ill , 

FST 

10 

2 

Mil  101 .HZI101  ,H3ll0l,HktlSl,GKtl01,Gmil  .(1(111  ,«t  101  ,K3t  10  1 . 

FST 

2 

AltlOt.AZtlOt  .A  31101,  81(10  t,  BZ  I1SI.BS  till,  OANPUI 131  .OAHPSfll.  Ill, 

FST 

2 

OUCt5.tl.SSIGLCllO,13l.NNU,ILC.ISC.NEtAI,S5l. 

FST 

2 

S3  IGVLC  (13,131  ,SSIGVSC(10,13I  , 

FST 

2 

SURIfcO,35I.N  tkO)  .NH.COSK  SSI,  CONI,  SS1GSC 1 10, 1 !1 , 

FST 

IS 

2 

NFS  (31  , 01SaL  B  «GH  ,  3N0TLC  <  10,  131  , 

FST 

2 

BN OTSC 110*131, BVELSC (10. 1 3> ,  SVELLC (10,131 «BSTOPtll ,  SVEINAX  <1 1 

FST 

2 

,NSAT , I TEST (18,131 

FST 

OINENSION  SNOUR 1 331, SSUIkO.331 .FINN (k0«3St.FIkVlkl.3Sl,8NIURt 3*1, 

FST 

2 

BXN0UR(35>,SSV  <k 0, 35) .SVNOUR  (351 

FST 

20 

DATA  RHO/1.44/  ,  PI/3.1kl59Z6/ 

FST 

• 

FST 

•initialize. 

FST 

• 

FST 

N*IS 

FST 

29 

nmntH) 

FST 

AA9EA*AREA(NI 

FST 

RR*R(Nt 

FST 

OO5L0RF*OCLDBFS(N) 

FST 

hho«hoini 

FST 

30 

HH1*H1 INI 

FST 

HH2*H2(N) 

FST 

NH3*H3(N) 

FST 

HHA*HA(N> 

FST 

GGK*G«(N) 

FST 

35 

GGV*GV(N) 

FST 

KK|*K1«NI 

FST 

KK2*K2tN» 

FST 

KK3>«S(N) 

FST 

AA1* At  INI 

FST 

AO 

AA2*A2(N) 

FST 

AA3*A3(N) 

FST 

881*91  INI 

FST 

862*82  INI 

FST 

883*83 (Nt 

FST 

AS 

C0M3*RH0*VrS(IV)«VF$(IV)*NN*AAREA*RR/(0XSPL6*6N) 

FST 

SSN  *  1.0 

FST 

ISO*  *  0 

FST 

75 

OANPStN.IiU)  *  OUCUV.l) 

FST 

IF  (ITERATE  «EQ.  0)  GO  TO  90 

FST 

50 

NTRY  *  0 

FST 

YP  *  0.0 

FST 

AO 

NTRY  *  NTRY  ♦  1 

FST 

X  *  YP 

FST 

DAMPS(N,IiUI  ■  OUCUV.1>*1rA**OUC(IV»2»*X*M.772  ♦  1*M»DUC(  XVt 3 

FST 

55 

Z1»X  »  fc.00»3UC(IV.kl*X»«Z.0  ♦  9,kO*OUC«IV,S»»X»»S.O  » 

FST 

2 

Zk.O»OUCtI¥,61»X»»k.| 

FST 

90 

PHIS  •  SSIGVLC  (M.IHU)  «  SSIS VSCtN, IHU)  ■  B, 

FST 

i 


l 


49 


SUBROUTINE  FINSTAB 


74/7%  OPT*0  ROUNO-*/  TRACE 


FTN  4.6*460 


11/26/79  11. 25.10 


SHOT SC IN.INUI *9VELSC (N.  IHU)*8NOTLC I N. I1UI*BVELLC(N. INU)  «  0. 

FST 

59 

• 

FST 

60 

60 

•BEGIN  LOOP  OVER  SPREADING  ANGLE  (NUI. 

FST 

61 

• 

FST 

62 

IF  (ISC  .E3.  01  NNU  «  1 

FST 

63 

DO  400  INU  *  1 .NNU 

FST 

64 

• 

FST 

65 

6$ 

•BEGIN  LOOP  OVER  NAVE  FREQUENCE. 

FST 

66 

• 

FST 

67 

DO  200  IM»t.NM 

FST 

60 

MME*NE(IW,INU) 

FST 

69 

NESQ«NN£*NME 

FST 

70 

70 

TUNF  *  MEdV.lNUI/NPHI 

FST 

71 

BA  *  2. •DAVPU f IHU) #TUNF 

FST 

72 

A  *  I.  -  TUNF*TUNF 

FST 

73 

CA  «  SORT (A* A  ♦  BA*BA> 

FST 

74 

* 

FST 

75 

75 

•CONFUTE  EFFECTIVE  LIFT  CURVE  SLOPE. 

FST 

76 

• 

FST 

77 

OCLDBE  *  H*0  ♦  MH1*MME  ♦  MH2*ME$Q  ♦  HM3*N£SQ9MME  ♦  HHL*ME$Q*ME$Q 

FST 

70 

OCLOBE  «  D2L  OBE •ODCLDBF* 100./ 3. 1415920 

FST 

79 

• 

FST 

•0 

•0 

•CONFUTE  ANPtITUOE  OF  FIN  ANGLE  TO  STABILX2E0  ROLL. 

FST 

01 

• 

FST 

02 

KR  »  KK1  -  MESQ*KK3 

FST 

03 

KI  »  MME*K<2 

FST 

04 

BR  •  BB1  -  NESQ*BB3 

FST 

05 

05 

BI  *  MME*B 52 

FST 

06 

C0N4»KR*KR  *  KX*KX 

FST 

07 

CONS  «  BR*3R  ♦  BI*BI 

FST 

00 

BAONS  «  6G<*GGV*S0RT (C0N4/C0N5I 

FST 

09 

• 

FST 

90 

90 

AR  ■  AA1  -  ME SQ*AA3 

FST 

91 

AI  *  NME*AA2 

FST 

92 

• 

FST 

93 

BS  «  2.*OA9PSIN.INUI*TUNF 

FST 

94 

CS  *  SORT! A* A  ♦  BS*BS) 

FST 

95 

95 

CON 6  «  AR*BR  -  AX*U 

FST 

96 

CON 7  «  ARMI  ♦  AI*BR 

FST 

97 

CONS  *  AR»AR  ♦  AIH1 

FST 

90 

SAONCSP  *  SSN*CON3*OCLOBE*BAONS/ (CS*SQRT (CONO)I 

FST 

99 

• 

FST 

100 

100 

COSMSI  ■  (KR*C0N6  ♦  KX*CON7)  /SORT  ( C0N4*C0N0*  C9N5I 

FST 

101 

SINKSI  •  UI*C0N6  -  KR*C0N7»/S0RT(C0N4*C0N0»C0N5> 

FST 

102 

• 

FST 

103 

•CONFUTE  ROLL  REDUCTION  FACTOR  FOR  THIS  FREQUENCY. 

FST 

104 

* 

FST 

105 

105 

SONU  «  CA/CS/SQRT ( 1.  *2.*S  AONCSPM  t A*COSKSI* 

FST 

106 

2  BS* SINRSII/CS ) *SAONC5P*S AONCSPI 

FST 

107 

• 

FST 

100 

•CONFUTE  STABILIZED  ROLL  *  FIN  NOTION*  AND  FIN  VELOCITY  SPECTRA. 

FST 

109 

• 

FST 

110 

110 

SSU (IM.INUI  »  SUR(IN.INU)  •  SONU  •  SONU 

FST 

111 

SSVI IN.INJt  *  SSU (IN. INU > •ME SQ 

FST 

112 

FINFIIN.INUI  *  SSU (IN. INU)  •  BAONS  •  BAONS  /  CONS 

FST 

113 

FXNV(IH.INU)  «  FINHIINtINUI  •  NESQ 

FST 

114 

* 

FST 

115 

V 1 


SUBROUTINE  FINSTAB  7-/7*  OPT«0  ROUNO**/  TRACE 


FTN  *.*»4t0 


«€N0  OF  LOOP  OVER  NAVE  FREQUENCY* 


200  CONTINUE 


•DETERMINE  RMS  ROLL,  FIN  NOTION,  AND  FIN  VELOCITY  VALUES, 


CALL  ALGRNS  I NN«M,SSU(i«INUIt  SHOUR I INUI I 
CALL  ALGRNG  ( NN,M.SSV II. INUI , SVNOUR (INUI I 
CALL  ALGRNG  I NN.N,FINNI1. INUI  ,BN|UR IINUII 
CALL  ALGRNS  « NM vM,FINVIl , XNUt * BVNOUR I INUI I 


•STORE  L0NGCRESTE3  VALUES  (NO  ITERATION  OVER  ROLL  OANPING  IS  OONEI • 


IF  (ISC  .NE.  0)  GO  TO  210 
SSIGLCIN.INUI  ■  SORT ISN0O* IINUII 
SSIGVLC  IN,  INUI*  SORT  (SVNOUR  (INUI  I 
BMOTLC (N, INUI  ■  SORT ( SNOUR IINUII 
BVELLCOWI1UI  *  SORT  I OVNO  UR  I  INUI  I 
IF  (ITERATE  • EQ •  01  GO  TO  A00 
SSGVMLC  •  SS I GVLC I N, INUI  /NPHl 
PHIN  *  DAMPS ( N, INUI 

CALL  ITREQI SSGVWLCI ,  RETURNS! 6 0,40 0,90 01 


•BEGIN  SUNNING  OF  SHORTCRESTED  RESPONSE  OATA* 


210  PHIS  *  PHIS  ♦  COS1 ( INUI *SN OUR (INUI 

SSI GVSC IN* I NUI  «  SSIGVSC IN, I NUI  ♦  COSt I INUI •SVNOUR ( INUI 
BNOTSC(N,INUI  «  BNOTSCIN, INUI  ♦  COS II INUI «SN0URI INUI 
BVELSC(N.IIU)  «  BVELSCIN* INUI  ♦  COS II INUI •SVNOUR I INUI 


SEND  OF  LOOP  OVER  NU. 


400  CONTINUE 

IF  (ISC  .El.  01  GO  TO  060 
PHIS  «  SORT (CONI* PHIS I 
SSIGVSC IN, I NUI  ■  SORT(CONl*SSXCVSCIN,IHUII 
8N0TSC (NfINUI  ■  SORT (CONI* BNOTSCIN, INUI I 
3VE LSC I N, IiUI  «  SORT (CONl*BVELSCIN, INUI I 
IF  (ITERATE  «EO*  01  GO  TO  250 
SSGVOM  *  SSIGVSC IN, INUI /NPNI 
PHIN  «  DAMPS (N, INUI 

CALL  ITRE0I  SSG VOMI ,  RETURNS! 80,250,9001 
250  SSIGSC IN, INJI  •  PNIS 

IF  (INSAT.EQ.OI  .OR,  (ISOR.EQ.il  I  GO  TO  000 
XI  *  (9ST0PI I VI /BVELNAXIIVI !•( SVELSC IN, INUI /BNOTSCIN, INUI I 
X2  «  SIN ( XII 
F  «  (4/PII»(X2/Xl) 

V  »  BSTOPIXVI  t  IBNOTSC (N, INUI *1 . 414211 

SSN  *  1/It  -  X2I  Mil  -  F*X2I  •ERRFIYI  ♦  IF  -  1 1 *X?*ERRF( Y/X2I  I 
IF  ISSN  .961  GO  TO  078 

ISOK  *  1 


GO  TO  75 

67  0  ITESTIN.IHJ! 
000  RETURN  AAA 
900  RETURN  BOB 
END 


FUNCTION  EIRE 


74/74  OPT*0  ROUNO* •/  TRACE 


FTN  k.  inti 


11/M/M 


1  c  erf 

REAL  FUNCTION  ERRFIX1  ERF 

oata  p/.«to  a//  'Ci/.staa?/  .cz/*.o«sti/  .ctmuh/  erf 

IF  (ASSIXI  .GT.  IS. I  GO  TO  5  ERF 

5  T  •  1/fl  ♦  »»XI  ERF 

C  ERF 

ERRF  «  1  -  (C 1*T  ♦  C2*T*T  ♦  CJ»T»T»T»*EXP«-*M»  ERF 

C  ERF 

GO  TO  10  ERF 

10  S  FRRF  =  1.  ERF 

10  RETURN  ERF 

ENO  ERF 


APPENDIX  C 
SPECIAL  ALGORITHMS 


ITERATION  OVER  ROLL-ROLL  DAMPING  EQUATIONS* 

If  roll  damping  coefficient  n  is  independent  of  roll  angle  (as,  for 
example,  in  some  cases  where  bilge  keels  are  not  appended  to  the  hull), 
FINCON  executes  with  a  unique  value  of  n  for  each  required  ship  speed. 
However,  as  has  often  been  found  through  model  experiments,  roll  damping  is 
dependent  on  roll  angle,  as  well  as  ship  speed  and  natural  roll  frequency. 
Thus,  the  program  requires  a  different  input  description  of  roll  damping. 

As  described  in  the  section  of  this  report  entitled  Program  Input,  n  is, 
in  this  case,  defined  by 

n1R  -  dQ  +  1.61  dqy0,772  +  l.M^y  +  4d,,y2  +  9.4d3y3  +  24d4y4  (7) 


where  y  ■  O'/ui^  and  the  fractional  power  q  *  0.772  arises  from  the  tur¬ 
bulent  skin  friction  contribution. 

Thus,  the  computational  problem  is  then  to  solve  the  roll  equation  of 
motion  involving  <J>  for  the  correct  value  of  nIR.**  This  has  been  accom¬ 
plished  by  finding  the  intersection  of  the  known  curve  of  Equation  (7) 
(i.e.,  n^R  *  f  (y)  where  y  =  a^ui(f))  and  the  initially  unknown  curve  g  * 
g(nIR>’  *s  computed  from  the  solution  of  the  roll-rate  equation  of 

motion,  using  the  roll  damping  value  nTD. 

The  solution  is  found,  in  brief,  by  the  following  procedure  as  applied 
to  the  example  of  Figure  4: 

1.  Assume  y  *  0,  determine  n  =  ffO) 

o  o 

2.  Using  nQ,  solve  the  roll-rate  equation  of  motion  for 


*o 


Of 

_£ 


g(nc) 


♦Taken  from  future  report  already  under  preparation  by  Cox. 

♦♦The  subscript  IR  is  used  here  since  it  is  assumed  that  motion  is 
taking  place  in  irregular  waves. 


n  -  «Iq  +  cLy**  +  dj  y  +  d2V2  +  dgy®  +  d^y* 


q  -  0.772 


ROLL  EQUATION 
OF  MOTION 


vi 


-C- 


np+1  np  "p-1 

ROLL  DECAY  COEFFICIENT,  n 

Figure  4  -  Illustration  of  Iteration  over  Nonlinear 
Roll-Roll  Damping 


3.  Assume  *  gQ/2,  and  determine  n^  *  f(y^) 

4.  Compute  g^^  *  gO^) 

5.  For  p  *  2,  3,  4  ...  etc. 

Calculate  yp  -  (gp_2yp_rgp.1yp.2)/[  “  (gP-r8P-2)] 

and  determine  n  ■  f(y  ) 

P  P 

6.  Compute  gp  *  g(np)  and  test  if  |l  -  gp/yp |  e  where  e  is  some 
small  number.  If  the  test  is  not  satisfied,  proceed  to  the  next 
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step  of  the  iteration.  If  the  test  is  satisfied,  then  the  iteration 
is  terminated  and  n^,  g  is  the  required  solution. 

In  FINCON,  e  is  usually  taken  as  0.01  and  is  called  PRCN  in  the  pro¬ 
gram  itself.  This  means  convergence  is  attained  when  |l  -  g^/y^l  £  0.01. 
Table  6  provides  a  typical  printout  of  the  steps  when  1PRINT(2)  is  input  as 
ITERATN.  IV  Indicates  the  speed,  such  as  the  first;  IMU  indicates  the 
heading,  such  as  the  seventh  or  90  degrees;  NTRY  indicates  the  numbers  of 
attempts;  PHIN  indicates  the  latest  selected  n  value;  YP+1  the  latest 
calculated  roll-rate  value  (divided  by  U),);  and  g  ,  y  ,  g  .  and  y  » 

<p  p— 1,  p— i ,  p—z,  p — / 

are  the  previous  calculated  roll-rate  values  (divided  by  u >^) ,  going  back 
in  time.  It  should  be  noted  that  the  printing  occurs  at  various  steps 
within  each  repetition  of  the  iteration  rather  than  just  at  the  end  (or 
beginning)  of  each  attempt.  In  the  particular  case  listed  in  Table  6,  the 
first  seven  lines  refer  to  the  unstabilized  case,  while  the  last  five 
refer  to  the  stabilized  case.  Also,  it  is  recalled  that  the  iteration  is 
operating  over  either  the  long-crested  or  short-crested  RMS  roll  rate. 

A  similar  procedure,  implemented  in  the  Navy  Standard  Ship  Motion  Program 
(SMP-79)  (six-degrees-of-freedom)  operates  over  the  resonant  region  of  the 
(long-crested)  roll  transfer  function.  Also,  the  roll-rate  values  actually 
used  here  for  the  internal  testing  (e.g.,  in  step  5  above)  are  usually 
taken  to  be  the  statistic  which  corresponds  to  the  experimental  data 
described  by  Equation  (7).  In  all  cases,  the  RMS,  single-amplitude  roll 
rate  is  used.  The  second  page  of  Table  4  provides  a  further  illustration 
of  this  example. 

GENERALIZATION  OF  COSINE  SQUARED  LAW  FOR 
SHORT-CRESTED  SEAS 

In  general,  a  15-degree  cosine  square  spreading  function  about  +  90 
degrees  is  used  for  calculation  of  ship  motions  in  short-crested  or  multi¬ 
directional  seas  (e.g.,  see  Reference  3).  However,  it  has  been  found  in 
calculations  done  previously,  that  a  more  refined  angular  spread  may  be  re¬ 
quired  for  roll-motion  calculations  at  higher  ship  speeds.  This  is  due  to 
the  highly  tuned  nature  of  roll  motion.  That  is,  when  considering  RMS 


! 

I 


! 

v 


Ik 


roll  motion  across  all  ship  he  *'  —  at  high  speeds,  the  maximum  values 
may  occur  in  between  the  standard  15-degree  increments  of  heading;  and,  if 
not  considered,  some  loss  of  energy  may  be  noticed. 

Therefore,  a  FINCON  program  option  is  available  for  varying  the 
spreading  angle  used.  Spreading  angles  of  5,  10,  or  15  degrees  may  be 
specified  on  data  card  9.  Though  the  spreading  is  thus  varied,  the  pre¬ 
dominant  heading  angles  are  never  varied  from  the  usual  15-degree 
increments.  Table  4  presents  results  when  a  15-degree  spreading  is 
specified.  In  general,  either  a  15-degree  spreading  (e.g.,  for  lower 
speeds)  or  a  5-degree  spreading  (e.g.,  for  higher  speeds)  is  recommended; 
although  no  specific  guidelines  for  the  use  of  either  are  currently 
available.  A  major  difference  between  the  two  angles  is,  of  course,  in 
program  run  time  and,  thus,  in  cost.  A  typical  15-degree  spreading  run 
may  cost  as  much  as  68  percent  less  than  a  5-degree  run.  In  general,  10- 
degree  spreadings  are  not  recommended,  because  irregular  trends  across 
ship  headings  may  be  perceived.  For  example,  for  a  75-degree  predominant 
direction,  the  10-degree  spreading  angles  are  at  -15,  -5,  5,  .  .  . ,  155, 

165  degrees;  and,  for  a  90-degree  predominant  heading,  they  are  at  0,  10, 
20,  .  .  .,  165,  180  degrees;  while  for  105  degrees  they  are  at  15,  25, 

35,  ....  185,  195  degrees.  If  roll  is  highly  tuned,  then  adjacent  short- 
crested  RMS  roll  angles  may  appear  to  have  erratic  behavior  (e.g.,  decrease 
at  75  degrees.  Increase  at  90  degrees,  decrease  at  105  degrees,  etc.). 

This  is  due,  numerically,  just  to  the  difference  in  the  base  spreading 
angles  in  each  case. 

The  algorithm  implemented  in  FINCON  to  provide  this  generalization  in 
the  short-cresting  procedure  is  defined  by  the  following: 

1.  Let  I  be  the  specified  spreading  angle  (5,  10,  or  15  degrees) 
or  angle  of  constant  energy 

2.  Set  i  -  (180/1) /2 
and  c  •  1/f. 

3.  Then, 

°sc  <W>  "  c  £  aLn(W  +  F)  cos2  F 

p— (Jl-1) 
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where  o£c  ■  squared  long-crested  RMS  roll  or  variance 
y  *  predominant  heading  angle 
2 

a  ■  squared  short-crested  RMS  roll  or  variance, 
s  c 

The  only  Interval  required  by  the  algorithm  is  1,  the  angular  interval  over 
which  a  constant  wave  energy-ship  response  is  assumed. 
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